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ABSTRACT 


The  design  of  piezoelectric  transducers  often  involves  the  satis¬ 
faction  of  multiple,  conflicting  design  specifications.  The  use  of 
conventional  mathematical  programming  techniques  in  this  endeavor  has 
been  inhibited  by  their  inability  to  consider  these  multiple  objec¬ 
tives  within  their  design  formulations.  Goal  Programming  is  an 
effective  alternative  methodology.  Transducer  design  analysis  through 
goal  programming  provides  the  optimal  solution  to  the  multiple  cri¬ 
teria  of  various  types  of  transducers.  The  solutions  to  five  design 
specifications  for  a  reasonably  complex  transducer  are  presented.  The 
details  of  the  formulation  and  analysis  of  results  provide  an  indica¬ 
tion  of  the  ability  of  this  approach  to  resolve  even  more  complex 
transducer  design  problems. 
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CHAPTER  1 


INTRODUCTION 


1. 1  Purpose 


This  study  provides,  through  Goal  Programming,  a  computer  tech¬ 
nique  to  generate  suitable  design  parameters  for  piezoelectric  trans¬ 
ducers,  given  a  set  of  objective  criteria.  The  specific  objectives  of 
this  thesis  are: 

1.  To  develop  a  mathematical  model  which  describes  the  trans¬ 
ducer  and  related  environment  and  which  accurately  predicts 
a  number  of  specific  responses  of  the  system. 

2.  To  adapt  this  model  to  the  Goal  Programming  formulation. 

3.  To  investigate  the  behavior  of  Goal  Programming  applied  to 
examples  of  "Tonpilz"  design. 

4.  To  provide  a  computer  program  to  enable  and  encourage  the  use 
of  the  technique  in  the  resolution  of  complex  transducer 
design  problems. 

1. 2  Problem  Definition 

One  style  of  transducer  favored  for  high-power  underwater-sound 
generation  is  the  "Tonpilz,"  a  composition  of  two  German  words — sound 
and  mushroom,  so  named  for  its  characteristic  shape  as  shown  in  Fig¬ 
ure  1.  Its  particular  construction  is  adaptable  to  the  requirements  of 
mounting  and  vibration  isolation.  The  design  of  these  transducers  involves 
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typically,  head  and  tail  masses  which  are  used  for  impedance  matching 
and  a  central  piezoceramic  motor  section.  The  actual  construction  of 
this  basic  three-section  model  sometimes  requires  additional  parts, 
such  as  insulators  and  electrodes,  as  well  as  an  assembly  lip  on  the 
head  called  the  "shank."  Thus,  a  typical  transducer  may  be  comprised 
of  different  materials  bonded  together  as  illustrated  in  Figure  1. 

Plane  longitudinal  stress  waves  passing  through  this  conglomerate 
travel  with  different  speeds  within  the  different  materials  and  suffer 
reflections  at  each  interface.  The  properties  of  these  materials 
must  be  carefully  considered  in  the  design  to  achieve  a  desired  re¬ 
sponse.  The  length,  weight,  resonance  frequency,  mechanical  quality 
factor,  impedance,  or  node  position  may  be  especially  important  in  a 
particular  application.  However,  satisfaction  of  one  requirement  may 
adversely  affect  another.  For  example,  the  design  of  a  transducer 
resonant  at  a  specific  frequency  and  of  minimum  weight  may  result  in  an 
unacceptable  length  or  an  undesirable  quality  factor.  Thus,  the 
design  of  a  transducer  for  a  specific  application  becomes  a  problem  of 
finding  an  optimal  compromise. 

One  previous  approach  has  utilized  mathematical  models  to  predict 
the  quality  factor  of  and  node  position  within  a  transducer  at  given 
frequencies  and  physical  dimensions.  The  consideration  of  other 
properties  has  not  been  possible  due  to  the  complexity  of  their  rela¬ 
tions  with  one  another.  The  need  for  a  flexible,  efficient,  multi¬ 
objective  optimization  technique  is  clear.  This  need  is  admirably 
answered  through  Nonlinear  Goal  Programming  (NLGP) ,  as  developed  by 
James  P.  Ignizio  [1].  The  specific  nonlinear  method  used  is  a 
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modification  by  Ignizio  of  a  well-known  method  of  nonlinear  single¬ 
objective  programming  known  as  Pattern  Search  [1] . 
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E  C  T  T 

s33*  ^33’  ^33*  an<*  £33  Pertain*  F°r  other  configurations  or  other 
types  of  crystals,  numbers  denoting  other  orientations  may  be  in 
order  [2,  9]. 

A  method  by  which  the  system  may  be  analyzed  begins  with  the 
equivalent  circuit  representation  for  the  ceramic  known  as  "Mason's 
Equivalent  Circuit."  This  circuit  displayed  in  Figure  2  shows  clearly 
two  mechanical  ports  and  one  electrical  port  [3,  4].  The  symbols  used 
in  Figure  2,  as  well  as  those  in  all  subsequent  figures  and  equa¬ 
tions,  are  defined  in  the  List  of  Symbols. 

With  W  defined  as  a  geometrically  dependent  function  of  length, 
the  circuit  elements  of  impedance  Z,  capacitance  Cq,  and  transformer 
turns  ratio  <j>  can  be  represented  as : 


(2) 

(3) 

(4) 

(5) 


This  circuit  can  represent  a  stack  of  individual  ceramic  segments  if 
the  | k5l |  of  each  segment  is  very  small  compared  to  it. 

To  represent  all  the  remaining  sections  of  the  transducer, 
multiple  T-sections  are  connected  smoothly  in  cascade  fashion  to  the 
ceramic  representation  at  the  mechanical  ports.  The  T-sections  for 
inert  materials  as  shown  in  Figure  3  are  equivalent  to  wave  guides. 


igure  2.  Mason's  Equivalent  Circuit  for  Piezoelectric  Material. 


Section  Circuit  Representation  for  Inert  Material 


.th- 


Equations  for  the  impedances  Z^  and  Z 2  of  the  i  T-section  are  given 
by: 


and 


Z1  =  ^PiCiAi 


Z2  =  -jP1ciAi/sin(kii.i) 


(6) 


(7) 


Note  the  identical  forms  of  Equations  (2)  and  (6)  and  of  Equations 
(3)  and  (7). 

On  placing  a  load,  Z^,  at  one  end  of  the  T-section  as  shown  in 
Figure  4,  the  input  impedance  as  seen  from  the  open  terminals  can  be 
calculated  using  Equation  (8)  : 

^ i t Zt  +  tan(k  l  )] 

in  [C±  +  jZL  tan(ki£i)]  *  (8) 

where  The  derivation  of  Equation  (8)  is  given  in 

Appendix  A. 

Use  of  the  wave-guide  equation  makes  it  particularly  easy  to 
remove  a  section  of  the  transducer  if  desired.  By  setting  the  length 
of  an  unwanted  section  equal  to  zero,  the  load  is  passed  unaltered  to 
the  next  section  as  can  be  readily  seen  by  putting  =  0  in  Equation 
(8),  viz.,  Z^  =  Z .  Figure  5  depicts  a  typical  transducer  represented 
by  Mason's  Equivalent  Circuit  modified  to  include  the  nonceramic 
sections. 

As  a  practical  matter,  the  transducer  must  be  mounted  in  a 


housing  and  also  coupled  to  the  acoustic  medium  (water,  in  this  case). 
A  common  practice  is  to  press  the  face  of  the  head  against  an 


Impedance  of  a  T-Section  with  Load  Z. 


Figure  5.  Transducer-Equivalent  Circuit  with  Mechanical  and  Electrical  Ports. 


acoustically  transparent  membrane  of  rubber  or  urethane  which  separ¬ 
ates  the  head  from  the  water.  The  pressure  is  applied  through  a  plate 
(commonly  called  a  web)  which  slips  down  over  the  body  of  the  trans¬ 
ducer  and  bolts  to  the  housing.  Figure  6  represents  a  transducer 
mounted  in  a  housing.  The  tail  is  generally  kept  free  to  vibrate 
into  air.  For  a  head-mounted  transducer,  a  compliant  support  structure 
is  introduced  between  the  web  and  the  back  of  the  head.  This  is  re¬ 
ferred  to  in  Figure  6  as  the  pressure-relief  material. 

The  entire  assembly  consisting  of  the  transducer  sections,  rubber 
membrane,  pressure-relief  material,  web,  and  loading  effect  of  the 
medium  can  be  represented  by  the  circuit  shown  in  Figure  7.  This 
circuit  is  a  fundamental  model  of  a  "Tonpilz"  transducer  and  its 
associated  housing.  Any  variation  in  radiation  loading,  mounting, 
or  electrode  configuration  can  be  incorporated  by  simply  altering  this 
circuit  representation. 

2. 2  Circuit  Solution 

The  input  impedance  is  evaluated  by  means  of  Equation  (8) , 
the  wave-guide  equation,  beginning  from  the  outermost  T-sections 
inward  toward  the  ceramic  section.  The  circuit  can  then  be  simpli¬ 
fied  to  the  one  shown  in  Figure  8,  where  Z[*-]  and  Z[->-]  are  the 
resulting  impedances  of  all  previous  T-sections  to  the  left  and 
right  of  the  ceramic  sections,  respectively. 

As  is  evident  from  Figure  8,  the  mechanical  impedance  presented 
to  the  transformer  can  be  written  as: 


Transducer  Mounted  in  Its  Housing 


Figure  7.  Circuit  Representation  of  Transducer  in  Housing. 


and  Z[+]  and  Z[->-]  represent  left  and  right  loads. 
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Z  ^mechanical  '  Z2  +  ‘Z1  +  Zl*»  >Z1  +  !W1/8W  +  ZM  +  2ZJ  . 

(9) 

The  input  electrical  admittance^  is  the  sum  of  the  electrical 
admittances,  l/R^  and  j2iTfCQ,  and  the  reciprocal  of  the  mechanical 
impedance  reflected  through  the  transformer: 

V  ’  1/K0  +  J2"fC0  +  *2/Z'f>mcha»ical  •  <10> 

where  f  is  the  frequency  in  hertz. 

The  term  containing  1/Rq  can  be  neglected  because  R^  is  normally 
very  large.  Equation  (10)  is  the  solution  to  the  circuit  representa¬ 
tion  and  is  used  to  calculate  the  electrical  admittance  at  any 
frequency.  The  impedance  Z(f )mechanicai  a  function  of  frequency 

through  the  wave  number  k  of  each  section. 

2.3  Stress-Bolt  Consideration 

The  glassine  molecular  form  of  the  ceramic  shatters  easily  under 
tension,  making  it  imperative  that  the  ceramic  be  kept  in  a  state  of 
constant  compression,  particularly  while  subjected  to  high-electric 
fields.  To  accomplish  this  compression,  a  bolt  is  passed  through  the 
ceramic  center  from  shank  to  tail.  Often  included  between  the  tail  and 
the  bolt  head  are  springs  called  Belleville  washers.  The  stress  rod 
(neglecting  Belleville  washers)  can  be  modeled  in  Mason's  Equivalent 
Circuit  as  another  T-section  joined  to  the  shank  and  the  back  of  the 
tail  [7].  Figure  9  represents  a  drawing  of  this  circuit. 


Transducer  Circuit  with  Stress  Bolt 
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Evaluation  of  the  circuit  indicated  in  Figure  9  is  difficult  using 
the  cascade  Equation  (3)  because  the  bolt  couples  the  elements  together. 
Falling  back  on  standard  circuit  analysis,  the  admittance  is  determined 
and  is  found  to  have  a  direct  functional  dependence  on  the  and 
elements  in  each  T-section.  Details  of  the  admittance  equation  of 
this  circuit  are  given  in  Appendix  B. 

The  difficulty  with  this  admittance  formula  appears  in  the  terms 
involving  Z ^  =  -jZ^/ sin(k^fc^) .  The  ability  to  eliminate  any  section 
by  letting  its  length  equal  zero  is  lost;  for  here,  a  zero  length 
generates  a  division  by  zero.  Programming  the  admittance  equation 
which  includes  the  stress  bolt  actually  limits  the  flexibility  of 
choosing  or  eliminating  sections  in  the  design.  Furthermore,  prac¬ 
tice  has  shown  that  the  addition  of  the  stress  bolt  in  the  transducer 
has  a  minimal  effect  on  the  resonance  frequency  and  admittance. 

Brickman  confirms  this  observation.  In  his  analysis  of  transducers, 
preloaded  with  stress  rods  and  Belleville  washers,  he  concludes  [10]: 

The  effect  of  changes  in  either  preload  (stress  bolt) 
or  isolator  stiffness  is  likely  to  be  rather  small  since  the 
Kp/Kq  and  Kr/Kq  stiffness  ratios  are  both  considerably  less 
than  unity  in  any  practical  transducer  design,  where 
Kp  =  stiffness  of  the  preload  rod 

Kr  =  stiffness  of  isolator  element  or  Belleville  washer 
Kq  =  stiffness  of  drive  element 

For  these  reasons,  the  stress  rod  has  been  neglected  in  the 
transducer  model,  and  Equations  (9)  and  (10)  have  been  chosen  to  mathe¬ 
matically  represent  the  transducer. 
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2.4  The  Admittance  Loop 

When  Equation  (10)  is  plotted  in  the  complex  plane  with  frequency 
as  a  parameter,  the  admittance  loop  shown  in  Figure  10  is  obtained. 

Five  important  design  parameters  can  be  calculated  from  the 
information  contained  in  the  admittance  loop.  They  are  the  quality 
factor,  the  electromechanical  coupling  coefficient,  the  first  resonance 
frequency,  and  the  magnitude  and  phase  of  the  admittance  at  resonance 
[5,  7].  These  quantities  are  listed  in  Table  1  with  a  brief  descrip¬ 
tion  of  their  importance  to  transducer  design. 


2.5  The  Node  Position 


The  approximate  position  of  the  velocity  node  at  f  in  the  transducer 
can  be  determined  using  a  wave-guide  analysis.  If  the  ceramic  section 
were  operated  without  an  external  load  on  either  face  (i.e.,  in  a 
free-free  boundary  condition) ,  it  would  behave  as  a  half-wave 
resonator  with  maximum  motion  at  the  ends  and  a  velocity  node  in  the 
center.  In  this  configuration,  the  node  represents  a  clamped  or  fixed 
point  and  the  ceramic  is  effectively  separated  into  two  independent 
quarter-wave  resonators  with  fixed-free  boundary  conditions. 

If  the  other  sections  are  added  to  the  ceramic,  the  boundary 
conditions  on  these  quarter-wave  resonators  change  from  fixed-free 
to  fixed-loaded.  To  derive  the  equations  to  predict  the  nodal 
position,  the  transducer  must  be  represented  by  a  wave  guide  with 
interfaces  between  the  loads  and  the  ceramic  at  A  and  B,  as  shown  in 
Figure  11. 


DIRECTION  OF 


Figure  10.  Admittance  Loop  of  Transducer  Circuit  with  Important  Frequencies  Marked. 
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The  condition  of  continuity  of  force  and  velocity  demands 
equality  of  mechanical  impedances  at  the  interfaces.  The  impedance 
looking  into  the  ceramic  from  A  toward  the  node  is  -jpcA  cot(kc&p). 
Therefore,  from  the  head/ceramic  interface  at  A,  we  must  have: 

ZM  =  "jPcA  cot  (k^)  .  (11) 

Similarly,  it  must  be  that  at  B,  the  relation  is: 

ZM  =  -jpcA  cot  (kc&B)  ,  (12) 

where  k  =  wave  number  in  the  ceramic, 
c 

Equations  (11)  and  (12)  are  easily  solved  for  2_  and  2_  after 

D  r 

calculation  of  the  quantities  ZM  and  ZM  from  the  repeated  applica¬ 
tion  of  the  wave-guide  equation  to  the  circuit  model.  Of  course, 
if  ZM  =  ZM,  then  2^,  =  2.^.  But  for  unequal  loads,  including  real 
components  (phase  differences) ,  2.^  /  2,^  and  the  node  becomes  merely 
a  minimum  velocity  point. 

It  is  highly  desirable  to  place  the  node  inside  the  ceramic 
rather  than  in  the  head  or  tail.  This  is  because  the  ceramic  develops 
an  EMF  in  response  to  a  stress;  and  since  the  stress  is  proportional 
to  the  derivative  of  the  velocity,  the  greatest  stress  occurs  in  the 
vicinity  of  the  velocity  node. 

2.6  Summary  of  Model 

The  transducer  and  its  associated  housing  are  represented  by 
T-sections  cascaded  onto  Mason's  Equivalent  Circuit  (Figure  7). 


The  circuit  is  solved  for  its  electrical  input  admittance  by  means  of 
Equation  (8).  This  admittance  is  repeatedly  evaluated  for  increasing 
frequency;  and  as  the  admittance  loop  is  traced  out,  its  frequencies 
of  interest  are  selected.  Table  1  presents  the  attributes  of  the 
transducer  which  are  evaluated  from  these  frequencies  and  concomitant 
admittances.  The  position  of  the  velocity  node  is  determined  by  the 
solution  of  the  circuit  in  the  same  manner  inward  to  the  ceramic  and 
application  of  Equations  (11)  and  (12) . 

2. 7  A  Verification  of  the  Computer  Algorithm 

In  order  to  assess  the  validity  of  the  computer  algorithm  used 
to  predict  the  correct  resonance  frequency,  a  three-section  transducer 
was  postulated  as  a  test  case.  Sections  one  and  three  were  composed  of 
brass,  while  section  two  was  piezoceramic.  All  sections  had  equal 
cross-sectional  areas  and  were  of  equal  lengths.  Such  a  simplification 
could  represent  an  unmounted  transducer  in  vacuum.  To  test  the 
algorithm  for  frequency,  the  parameters  of  the  hypothetical  three- 
section  transducer  as  shown  in  Table  2  were  inserted  into  the  input . 

The  other  available  sections  of  the  design  model  were  removed  by 
letting  their  lengths  equal  zero.  A  small  loss  factor  of  0.001  was 
permitted  in  the  ceramic,  and  the  web  and  radiation  loads  were  set  to 
zero  corresponding  to  an  unmounted  transducer  with  air  (vacuum) 
loading.  The  magnitudes  of  the  clamped  capacitance,  CQ,  and  the  trans¬ 
former  turns  ratio  <|>  were  arbitrary,  as  these  do  not  affect  the  value  of 
the  frequency  at  the  fundamental  resonance.  Using  these  data,  the  computer 


algorithm  evaluated  the  admittance  and  selected  the  frequency  of 


maximum  conductance,  the  resonance  frequency,  as  f^  =  21770  Hz. 

It  can  also  be  shown  that  the  resonance  frequency  can  be  calcu¬ 
lated  by  demanding  equality  of  impedance  at  any  interface.  This 
leads  to  the  development  of  a  relation  for  any  three-section  composite 
half-wavelength  longitudinal  vibrator: 


tan(k2^2) 


tan  tan  (k^^) 

tan  (k.^)  tan  (k.^)  ~  1 


(13) 


where 


and 


P1C1A1 
B.  =  ■ 

1  P2°2A2 


=  P3C3A3 
3  P2C2A2 


with  A  =  area,  p  =  density,  c  =  bar  velocity,  k  =  wave  number,  and 
H  =  length. 

The  Secant  Method  was  chosen  to  find  the  roots  to  Equation  (13). 
From  the  data  shown  in  Table  2,  the  resonance  frequency  of  the  thtee- 
section  test  case  given  by  this  method  was  21769.221  Hz.  This  is  in 
excellent  agreement  with  the  value  obtained  from  the  computer  algorithm 
and  so  lends  confidence  to  the  use  of  this  routine. 

As  is  well  known  [9] ,  the  presence  of  losses  in  a  parallel 
circuit  decreases  the  resonance  frequency.  The  ceramic  transducer 
is  in  essence  a  parallel  circuit  as  shown  in  Figure  8  and,  therefore, 
should  suffer  a  decrease  in  its  resonance  frequency  for  increased 
losses.  To  display  this,  the  data  from  Table  2  were  rerun  with  the 


loss  factor  of  the  ceramic  section  set  equal  to  0.2  Instead  of 
the  previous  value  of  0.001.  The  result  Is  a  smaller  resonance 
frequency  as  expected: 


f  =  21560  Hz  . 

y 


2.8  A  Verification  of  the  Model 


Five  quite  different  transducers  constructed  with  a  variety  of 
materials  were  selected  to  serve  as  inputs  for  the  algorithm  which 
evaluates  the  transducer  model.  Physical  dimensions  of  each  trans¬ 
ducer  were  measured,  and  tests  using  the  Vector  Admittance  Locus 
Plotter  (VALP)  were  made  on  each  to  determine  the  frequencies  of 

resonance.  Quality  factors  Q  were  obtained  from  f,  and  f„,  being 

m  1  z 

the  frequencies  of  maximum  and  minimum  susceptance,  respectively,  and 
from  the  resonance  frequency  f  by: 


(14 


The  effective  coupling  coefficient  was  estimated  from  the  approxi 

mate  equation: 


(15 


where  f  is  the  approximate  frequency  of  minimum  conductance. 
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Values  of  the  magnitude  of  the  admittance  at  resonance  were  re¬ 
corded.  In  two  cases,  the  position  of  the  node  was  measured  by 
placing  an  accelerometer  probe  on  the  ceramic  and  searching  for  the 
position  of  minimum  longitudinal  movement. 

Table  3  shows,  for  the  five  differently  constructed  transducers, 
the  percentage  differences  between  the  values  returned  by  the  program 
and  those  measured  experimentally. 

’■  On  averaging  the  first  four  values  in  column  1,  it  can  be  seen 
that  the  program's  calculated  frequencies  are  approximately  7.1% 
higher  than  the  mean  measured  resonance  frequencies.  In  order  to  see 
whether  this  discrepancy  could  be  ascribed  to  such  neglected  factors 
as  the  finite  thickness  of  the  electrodes  and  glue  joints,  additional 
lengths  of  1.0  mm  were  added  in  the  program  to  each  pf  the  insulators. 
However,  this  adjustment  only  lowered  the  calculated  resonance  fre¬ 
quencies  by  an  average  of  0.3%.  Therefore,  it  is  thought  that  other 
mechanisms  such  as  transverse  vibrations,  head  flap,  or  compliance 
of  the  glue  must  account  for  this  shift.  The  disparate  value  of  a 
14.4%  difference  for  transducer  #5  was  assumed  to  be  caused  by  a  defect 
in  this  transducer.  It  is  not  uncommon  for  a  transducer  to  have  a 
poorly  glued  joint  which  manifests  itself  as  a  widely  deviated  reson¬ 
ance  frequency,  or  even  quite  often  as  two  resonances.  In  fact,  one 
of  the  tests  for  this  condition  is  to  compare  the  resonance  frequency 
with  a  sampling  of  others.  Insofar  as  this  aspect  is  concerned, 
transducer  // 5  may  presumably  be  discounted. 

Column  2  contains  the  percentage  differences  between  calculated 


and  measured  quality  factors.  Normally,  with  air  loading  on  both 
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Table  3 

Percentage  Differences  Between  the  Program-Generated 
Design  Values  and  the  Measured  Values 


Transducer  //  %Af 

y 

%AQ 

m 

%Akeff 

%A|  Admit  |  (m  ) 

%A£p 

1 

7.20 

-  6.47 

-57.00 

60.30 

-40.05 

2 

6.68 

31.60 

23.20 

-58.10 

-  6.10 

3 

8.71 

9.47 

-  3.53 

60.50 

NA 

4 

5.98 

-11.07 

-  8.03 

22.70 

NA 

5 

14.40 

40.10 

-56.00 

75.00 

NA 
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acoustic  ports,  a  transducer's  Q  will  lie  between  200  and  500.  The 

m 

differences  in  column  2  are  well  within  the  deviations  encountered  in  a 
set  of  "identical"  elements.  The  quality  factor  is  a  parameter  which 
is  greatly  dependent  on  the  real  part  of  the  external  load,  as  well  as 
on  the  internal  losses  of  the  system.  When  the  input  data  for  the  five 
transducers  were  rerun  using  water  loading  on  the  faces,  was  seen 
to  fall  to  values  ranging  between  5  and  9,  while  at  the  same  time,  the 
frequency  f^  shifted  very  slightly  downward.  This  behavior  is  con¬ 
firmed  by  experiment  [9]. 

The  coupling  coefficient  discrepancies  are  displayed  in  column  3. 
Here,  there  is  very  poor  agreement  between  calculations  and  measure¬ 
ment.  It  is  felt  that  the  measured  is  more  probably  in  error 

because  of  the  approximations  involved  in  the  measurement  technique 
using  the  VALP.  The  calculated  values  for  ke^  in  air  range  from 
0.50  to  0.65  and  are  probably  more  reliable.  In  addition,  for  water 
loading,  the  calculated  k  ^  increased  on  average  by  14%  in  agreement 
with  experimental  observations  [9] . 

Column  4  deals  with  the  absolute  value  of  the  admittance  at 
resonance.  These  percentage  differences  are  well  within  the  variability 
shown  by  a  group  of  "identical"  transducers,  and  agreement  is  considered 
to  be  good.  For  water  loading,  the  program-generated  values  for  the 
magnitude  of  the  admittance  in  all  cases  decreased  by  two  orders  of 


magnitude,  as  expected. 


The  last  column  of  Table  3  contains  results  for  the  nodal  posi¬ 
tion  in  two  of  the  transducers.  Agreement  was  reasonable  for  trans¬ 
ducer  //2,  but  rather  poor  for  #1.  Generally,  it  is  expected  that  the 
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calculated  position  of  the  node  will  be  accurate  enough  for  use  in  the 


design  program 


*+*&**+»* 


CHAPTER  3 

GOAL  PROGRAMMING 


3.1  General  Theory 

In  the  Nonlinear  Goal  Programming  formulation  (NLGP)  chosen  for 
this  design  problem,  the  independent  variables  are  assigned  to  the 
vector  of  unknowns  x.  The  dependent  variables  are  called  the  objective 
functions,  and  are  the  components  of  the  vector  Y.  The  initial  step  in 
the  goal-programming  formulation  is  the  heuristic  development  of  a  set  of 
mathematical  statements,  which  are  functions  of  the  independent 
variables  x  to  describe  the  various  design  objective  functions. 

i 

The  general  form  of  the  equation  relating  objective  function  to 
desired  value  is: 


Yi(x)  +  ni  -  p±  =  b±  .  (16) 

Here,  the  vector  element  b^  is  the  desired  achievement  or  optimum 
value  acceptable  for  Y^.  In  some  cases, the  number  b^  may  actually  be 
unattainable,  whereupon  a  reasonable  goal  must  be  established.  For 
example,  it  may  be  desirous  to  have  a  cost  objective  equal  to  zero; 
however,  a  more  reasonable  goal  of  20%  of  the  estimated  cost  should  be 
assumed . 

The  other  two  variables  introduced  in  Equation  (16)  are  the  de¬ 


sign  deviation  variables,  n^  and  p^.  These  nonnegative  variables  are. 
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respectively,  the  amount  by  which  is  under  or  over  its  goal,  b^. 

To  achieve  the  desired  right-hand  side,  certain  of  the  deviation 
variables  must  be  driven  to  zero.  To  strive  for  an  exact  fulfillment 
of  a  goal,  both  the  negative  and  positive  deviations  will  be  mini¬ 
mized,  whereas  if  the  goal  is  unattainable  only  one  of  the  deviation 
variables  will  be  necessary. 

All  of  the  pertinent  n^  and  p^  are  multiplied  by  weights  w^  and 
summed  in  priority  levels  to  create  the  achievement  function  a,  as 
defined  in  Equation  (17): 

),  .  .  .  ,  (17) 

- 

where  k  =  index  of  priority  level  1. 
j  =  index  of  priority  level  2. 

The  achievement  function  is  split,  vector  fashion,  into  priority 
levels.  The  deviation  variables  relating  to  the  primary  objectives, 
those  which  must  be  satisfied,  are  summed  in  priority  level  //l. 
Minimization  in  this  priority  level  is  tested  before  consideration  of 
the  next  lower  level,  and  this  level's  solution  can  never  be  degraded. 
The  remaining  n^  and  p^  from  the  secondary  functions  are  contained  in 
other  lower  priority  levels.  A  lexicographic  minimum  is  sought  for  the 
achievement  function.  Lexicographic  minimum  implies  an  ordered  mini¬ 
mum  in  which  the  priority  levels  are  considered  sequentially.  For 
example,  by  this  designation,  the  vector  (2,  20000)  is  a  smaller 
vector  than  (5,  1). 

Also  included  in  the  achievement  function  are  weight  factors. 

These  weight  factors  w  can  be  extremely  important  in  a  problem  with 


I  =  E  (wknk  +  wkPk) 
k 


E  (w.  n.  +  w.  p. 
•j  j  J  3*3 


35 


diverse  objectives.  Since  the  magnitude  of  n^  or  p^  may  differ  by 

several  orders  of  magnitude  for  different  objectives  at  the  same 

priority  level,  it  is  necessary  to  reflect  their  relative  importance 

in  the  achievement  function.  Accordingly,  a  weight  factor  is  assigned 

to  each  n^  and  p^  within  each  priority  level.  These  weight  factors 

also  provide  versatility  to  change  the  importance  of  one  or  another 

term  within  the  achievement  function. 

The  entire  NLGP  technique  for  optimization  using  two  priority 

levels  can  be  succinctly  expressed  as  follows  [1] : 

Locate  the  transducer  design  variables  x  so  as  to  lexico¬ 
graphically  minimize  a: 


+  V’k5*  f(Vj  +  wjPj> 

j  J 


such  that : 


Y^Cx)  +  n^  -  p^  =  b^  ,  for  all  i  =  1  .  .  .  m 

where:  k  =  summation  index  for  primary  objectives  in  priority 
level  //I , 

j  =  summation  index  for  secondary  objectives  in  priority 
level  //2, 

m  =  //  of  objective  functions  (m  =  k  +  j  ) . 

The  method  of  pattern  search  begins  by  changing  a  variable  by  a 
small  amount.  A,  and  evaluating  the  resultant  achievement  function.  If 
the  new  achievement  function  is  smaller  by  a  specified  tolerance  than 
the  old  one,  the  A  step  is  deemed  successful.  The  variable  is  then 
promoted  temporarily  to  the  new  value,  and  the  program  progresses  to 
test  the  next  variable.  At  the  end  of  a  successful  search,  all  vari¬ 
ables  are  accelerated  (each  in  the  direction  of  its  successful  A) 
by  an  acceleration  factor  a.  The  achievement  function  is  perturbed 
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at  this  accelerated  point,  and  if  the  result  is  smaller  than  before, 
the  process  begins  anew.  If  no  improvement  is  found,  the  search 
falls  back  to  the  previous  base  point,  the  size  of  each  A  is  reduced, 
and  the  search  is  begun  again. 

There  are  three  methods  for  reaching  convergence  to  an  optimum 
compromise.  The  first,  the  most  obvious,  is  to  completely  satisfy 
all  goals  and  reach  a  =  (0,  0).  Failing  that,  the  program  looks  at 
the  size  of  the  change  in  a  after  each  perturbation.  If  the  change 
is  less  than  0.001  in  each  priority  level,  then  the  program  assumes 
a  local  optimum  point  has  been  reached  and  declares  success.  If 
neither  of  these  states  occur,  the  program  stops  with  a  best-vector- 
to-date  for  x,  after  using  the  maximum  number  of  cycles  or  reductions 
the  operator  has  specified. 

The  NLGP  code  chosen  for  the  transducer  design  program  is  an  early, 
simple  version  of  Nonlinear  Goal  Programming.  More  powerful 
codes  exist  which  can  solve  problems  with  80,000  variables  and 
several  thousand  objective  functions  and  at  the  same  time  are  (typic¬ 
ally)  many  orders  of  magnitude  faster  while  requiring  less  storage  [11]. 
However,  for  this  small  problem,  the  simplified  version  is  adequate. 

3. 2  Application  to  Transducer  Design 

Applying  the  NLGP  technique  to  transducer  design,  the  sections 
of  the  transducer  and  housing  are  numbered  according  to  Figure  12. 

The  possible  independent  variables  x  are  (for  each  of  the  eight 
transducer  sections)  length,  cross-sectional  area,  density,  bar 
velocity,  and  loss  factor,  for  a  total  of  40  possible  variables. 


Transducer  Section  Numbers  Used  in  NLGP  Program  for  Design 
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Currently,  11  objective  functions  are  defined.  These  are  not 
the  only  objectives  which  might  be  chosen,  and  others  can  easily  be 
added  or  substituted.  The  goals  b^  are  chosen  as  typical  transducer 
specifications.  These  objective  functions  are: 


(x)  =  primary  admittance  resonance  frequency,  f  ,  and  b^ 


specification  value. 


Y2  (x)  =  ratio  between  the  lengths  of  the  ceramic  backward  and 

forward  of  the  velocity  node,  SL.  ,  /i.  ;  b„  =  1.0. 

back  front  2 

Y.,  (x)  =  mechanical  quality  factor,  Q  ;  b.  =  0.0. 
j  m  j 


Y^  (x)  =  electromechanical  coupling  coefficient,  k^^;  b^  =  1.0. 

6 

Y,.  (x)  =  total  transducer  length,  £  SL  ;  b,.  =  0.0. 

i=l  1  3 

6 

Y,  (x)  =  total  transducer  weight,  £  P.A.A.;  b,  =  0.0. 

o  1  i  1  0 

Y^  (x)  =  magnitude  of  the  admittance  at  resonance,  \u\  ^  ;  b^  = 

"  y 

0.1. 


Yg  (x)  =  phase  of  the  admittance  at  resonance, 
bg  =  0.0. 


Y  (x)  =  cost  of  materials,  £  A.£.  •  (unit)  .  bn  =  20%  of 

7  i  l  cost  1  9 

estimated  cost. 
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Yn(x)  =  ratio  of  the  length  of  the  ceramic  in  front  of  the  node 

to  the  total  length  of  ceramic,  /SL  ,  ;  h  = 

front  ceramic*  11 


These  objective  functions  are  evaluated  in  the  design  program  in  the 
subroutine  YVALUF  following  the  methods  and  equations  given  in 
Section  2.6. 

For  most  typical  "Tonpilz"  design  problems,  the  achievement  function 
is  written  as: 


3  ^W10P10  +  WllPll^wlnl  +  W1P1  +  w2n2  +  W2P2  +  W3P3 

+  w^n4  +  w^p5  +  w“p6  +  w^n7  +  WgPg  +  w“p9)  . 


This  formulation  shows  the  placement  of  the  deviation  for  objective 

functions  10  and  11  in  priority  level  #1.  Objective  function  Y^Cx)  is 

the  ratio  £,/£  .  ,  and  objective  function  Y, ,  (x)  is  £../£,  .  .  To 

f  ceramic*  J  11  b  ceramic 

insure  the  node  remains  inside  the  ceramic,  neither  of  these  ratios 
should  become  greater  than  unity.  Therefore,  the  b.^  of  these  objective 
functions  are  set  equal  to  unity.  There  is  a  definite  advantage  to 
placing  the  deviations  of  these  ratios  in  priority  level  one.  If  the 
starting  point  or  initial  guess  for  the  variables  is  chosen  so  as  to 
place  the  node  inside  the  ceramic,  then  the  nodal  ratios  are  less  than 
unity,  and  p^  and  p  will  both  be  zero  after  the  first  achievement 
function  test.  As  the  variables  are  perturbed  in  the  search  for  a 
lexicographic  minimum,  this  priority  level  solution  of  zero  cannot  be 
degraded,  so  that  now,  all  future  solutions  are  guaranteed  to  have  their 
node  inside  the  ceramic  boundaries. 


Table  4  displays  a  list  of  some  of  the  important  variables  of 
the  NLGP  code  for  transducer  design.  Table  5  provides  a  User's 
Guide  to  the  input  for  this  program.  The  entire  program  is  listed  in 
Appendix  C. 

3.3  Variables 

Any  of  the  parameters  such  as  density  P,  bar  velocity  c,  area  A, 
length  l,  or  loss  factor  ri  of  each  section  may  be  chosen  as  a  variable 
to  be  optimized  by  the  NLGP  technique.  For  proper  operation,  the 
selected  variables  must  be  placed  as  components  in  the  vector  x.  This 
takes  place  in  the  subroutine  "VARIB"  at  the  end  of  the  program  deck. 
For  example,  to  let  the  density  of  section  1  and  the  sound  velocity 
of  section  3  be  variables,  cards  containing  the  following  assignments 
must  be  placed  in  "VARIB." 

RH0(1)  =  X(l) 

C(3)  =  X(2) 

For  theoretical  research,  any  of  the  material  parameters  may  be 
variable;  however,  for  practical  design,  there  are  some  limitations. 

In  actual  practice,  the  product  of  density  and  bar  velocity,  the 
specific  acoustic  impedance,  is  the  quantity  of  interest,  and 
materials  exhibiting  the  necessary  combination  are  limited  by  avail¬ 
ability  and  cost  to  two  or  three  (in  each  section).  It  is  far  easier 
to  fix  the  specific  acoustic  impedances  throughout  and  then  to 
optimize  the  remaining  parameters  than  to  try  to  obtain  the  materials 
corresponding  to  a  program-generated  solution. 
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Table  4 

Definition  of  NLGP  Code  Variables 


Variables 


Utility 


NVAR 

NOBJ 

NPRIOR 

NACH 

NMAX 

X 

EPS 

RHS 

MPRI 

WEIGHT 

ALPHA 

BETA 

MAXRED 

IPRINT 

EPSY 

SIGN 


IROW 

Y 

XBASE 

ZTEST 

RHO 

C 

L 

A 

ETA 

TURNS 

CAP 

ZW 

ZPR 

STEP 


Number  of  Decision  Variables 
Number  of  Objective  Functions 
Number  of  Priority  Levels 
Number  of  Terms  in  Achievement  Function 
Maximum  Number  of  Pattern  Search  Cycles  Allowed 
Vector  of  Decision  Variables  at  Test  Point 
Vector  of  Step  Sizes  Used  in  Pattern  Search 
Right-Hand-Side  Values  of  Objective  Functions 
Priority  Level  of  Given  Deviation  Variable 
Weight  Factor  of  Given  Deviation  Variable 
Pattern  Search  Acceleration  Factor 
Step  Reduction  Factor 

Maximum  Number  of  Step  Reductions  Allowed 
Output  Switch — 0  Gives  Only  Final  Result; 

1  Gives  Result  of  Every  Pattern  Search  Cycle 
Vector  of  Improvement  Tolerances  in  Achievement 
Indicates  Positive  or  Negative  Deviation  to  be 
Minimized  from  Objective  Function  IROW  at 
MPRI  Level  with  Weight  Factor  Weight 
Objective  Function  of  Given  Deviation  Variable 
Vector  of  Functions  of  Decision  Variables 
Base  Point  for  Hooke-Jeeves  Algorithm 
Achievement  Function  Value  Vector  at  a  Test  Point 
Vector  of  Density  Values  (MKS)  KG/CU  A 
Vector  of  Sound  Speeds,  m/sec 
Length  of  Each  Wave-Guide  Section,  M 
Area  of  Each  Section  Square  Meters 
Vector  of  Loss  Factors 

The  Square  of  the  Transformation  Factor 
The  Calculated  Capacitance 
The  Loading  Presented  to  the  Front  Face 
The  Loading  Presented  by  the  Web  Applied  to  the 
Pressure  Relief 

The  Frequency  Increments  Used  to  Find  Resonance 
by  Iteration 

Vector  of  Values  of  Each  Cost  Per  Unit  Volume 


COST 
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Table  5 

A  User's  Guide  to  the  Input 


Variable 


Format 


Card  1 

NVAR,  NOBJ,  NPRIffR,  NACH,  NMAX  515 

Card  2 

IPRINT,  MAXRED  215 

Card(s)  3- 

X(I),  1=1,  NVAR  8E10.4 

Card(s)  4- 

EPS(I),  1=1,  NVAR  8E10.4 

Card(s)  5- 

RHS(I),  1=1,  NOBJ  8E10.4 

Card  6 

ALPHA,  BETA  2E10.4 

Card  7 

EPSY(I),  1=1,  NPRIOR  8E10.4 

Cards  8-22 

SIGN(I) ,  IRGW(I),  MPRI(I) ,  WEIGHT(I), 

1=1,  NACH  F5.0,  215,  F10.0 

Cards  23-24 

COST (I) ,  1=1,  8  8E10.4 

Cards  25-34 

RHO(I) ,  C(I),  L(I) ,  A(I),  ETA(I) ,  1=1,  8  8E10.4 

Card  35 

TURNS,  CAP,  ZW,  ZPR  6E10.4 

NOTE:  ZW  AND  ZPR  ARE  COMPLEX  AND  REQUIRE  2 

DATA  FIELDS  EACH 


Card  36 
STEP 


15 


While  at  first  glance  it  would  seem  advantageous  to  include  the 


cross-sectional  area  as  a  variable,  there  are  strong  objections  to 
this.  The  area  of  the  head  is  usually  fixed  by  array  dimension  and 
wavelength  requirements.  The  outside  diameters  of  the  shank,  insu¬ 
lators,  and  tail  are  usually  chosen  to  conform  to  the  outside  diameter 
of  the  ceramic  for  purposes  of  alignment.  In  turn,  the  ceramic  can  be 
purchased  most  inexpensively  if  it  is  selected  in  stock  sizes. 

Finally,  a.  most  important  consideration  is  the  evaluation  of  the 
capacitance,  C^,  which  strongly  affects  the  value  of  the  coupling 
coefficient  and  also  the  value  of  the  admittance  at  resonance.  The 
area  of  the  ceramic  presented  to  the  plane  wave  is  not  necessarily  the 
area  of  the  electrodes  which  define  the  capacitance.  There  are  many 
electrode  configurations  other  than  the  sandwich  type,  such  as  cylin¬ 
drical  shells  with  the  electrodes  inside  and  outside,  or  in  stripes 
down  the  sides.  Thus,  if  the  area  of  the  ceramic  face  is  to  be  a 
variable,  the  capacitance  must  be  made  to  assume  variability,  also, 
in  the  appropriate  manner.  No  attempt  has  been  made  to  program  all 
the  variety  of  capacitance  relations  possibly  applicable.  Instead, 
both  area  and  length  of  the  ceramic  should  be  fixed,  and  the  capaci¬ 
tance  calculated  as  required.  It  is  felt  that  this  method  will  pro¬ 
vide  the  most  flexibility  in  geometry  and  electrode  position  for  trans¬ 
ducer  design. 

Finally,  for  practical  purposes,  the  loss  factor  should  probably 
be  kept  constant  because  it  might  be  extremely  difficult  to  find 
materials  which  match  all  the  other  desired  properties  as  well  as  any 
specification  on  losses. 


A 
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3.4  Weight  Factors  in  the  Achievement  Function 

The  weight  factors  w^  (multipliers  of  the  deviation  variables  n_^ 
and  p^)  are  extremely  important  to  proper  program  operation.  In  order 
to  have  an  order  of  magnitude  equality  among  all  the  terms  of  the 
achievement  function,  it  is  necessary  to  examine  the  anticipated  size 
of  the  deviations.  The  goals  of  the  example  in  Section  4.1  are 
examined  in  Table  6  to  illustrate  the  assignment  of  weight  factors. 
Regarding  the  assignments  of  the  w^  shown,  there  are  two  points  which 
should  be  made.  First,  all  of  the  deviations  are  slowly  changing 
with  the  exception  of  the  one  associated  with  the  resonance  frequency. 
The  weights  assigned  for  this  deviation  will  emphasize  the  solution 
for  f  until  this  frequency  approaches  its  goals,  whereupon  the 
other  objectives  begin  to  assume  importance.  Secondly,  the 
second  objective  function,  ^jjack^front  >  is  t0  be  optimized  toward 
1.0,  and  failing  this,  it  is  preferred  that  the  ratio  be  smaller  than 
unity  so  that  the  node  will  be  more  toward  the  tail.  The  following 
argument  demonstrates  the  direction  of  the  emphasis  of  the  solution. 
When  the  fraction  is  less  than  1.0,  P2  is  zero  and  ^  can  only  take 
values  less  than  1.0.  If  the  ratio  is  greater  than  1.0,  n2  is  zero 
but  P2  can  assume  any  value.  Thus,  there  will  be  a  larger  influence 
on  the  summation  a  if  the  deviation  is  positive,  and  this  means  the 
solution  will  be  more  concerned  with  lowering  P2  than  raising  n2»  a 
desirable  result.  The  rest  of  the  weights  are  assigned  in  a  straight¬ 
forward  manner  according  to  the  anticipated  size  of  the  respective 


deviations . 


|  admittance  | C  0.01  ~  .009  1000 


Table  6  (Continued) 


CHAPTER  4 


EXAMPLE  OF  DESIGN  PROGRAM  USAGE 


4.1  An  Illustrative  Example 

To  display  the  inner  working  of  the  design  program,  the  follow¬ 
ing  example  was  chosen.  Suppose  one  wished  to  calculate  the  lengths 
of  the  head,  insulators,  shank,  and  tail  of  a  transducer  necessary  to 
provide  an  element  resonant  at  5000  Hz  in  water,  with  no  housing,  and 
with  a  combination  of  other  design  features  as  defined  in  Table  6. 

To  begin  the  problem,  decisions  concerning  the  compositions  of  the 
sections  and  the  length  and  area  of  the  ceramic  must  be  made.  The 
capacitance  and  transformer  turns  ratio  applicable  to  the  geometry 
and  material  constants  must  be  calculated  (see  Equations  (4)  and  (5)). 
Initial  lengths  for  the  various  sections  to  be  optimized  can  be 
assigned  by  estimating  the  total  transducer  length.  Assuming  it  to 
be  a  half-wavelength  vibrator  of  ceramic,  then: 

1  =  -fj  ,  (19) 

where,  f  =  the  desired  frequency  of  resonance, 
c  =  bar  velocity  in  the  ceramic. 

After  the  length  of  the  ceramic  (and  other  fixed-length  sections)  has 
been  subtracted,  the  remaining  length  can  be  divided  among  the  addi¬ 
tional  desired  parts. 


, 
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It  should  be  emphasized  that  the  purpose  of  this  study  is  to 
develop  and  illustrate  a  design  technique,  not  to  provide  specific 
transducer  designs.  For  this  example,  the  head  and  shank  were  chosen 
to  be  aluminum,  the  insulators  were  glass,  and  the  tail  was  tungsten. 

Chosen  for  the  ceramic  were  four  discs,  each  of  cross-sectional  area 

-3  -2 

equal  to  2.38  x  10  square  meters,  and  length  5.0  x  10  m,  and 

f  R\ 

having  the  material  constants  of  Channelite  #5400 v  '  [10].  The 
capacitance  was  calculated  to  be  1134  picofarads  and  the  square  of 

the  turns  ratio  was  0.5576.  The  initial  lengths  of  the  head  and  tail 

-2  -3 

were  5  x  10  m.  One  insulator  length  was  2.5  x  10  m  while  the 

-3  -3 

other  was  3.0  x  10  m.  The  shank  length  was  chosen  to  be  3.0  x  10  m. 

The  assignment  of  the  variables  in  the  subroutine  "VARIB"  was 
made  as  shown  in  Table  7. 

Table  8  contains  a  listing  of  all  the  input  variables  for  this 
problem.  Explanations  of  all  variables  can  be  found  in  Table  4. 

The  radiation  loading  was  set  to  the  value  of  the  specific  acoustic 
impedance  (pc)  of  water.  Since  housing  effects  were  to  be  neglected, 
the  web  load  and  the  lengths  of  sections  7  and  8  were  set  equal  to 
zero.  Cross-sectional  areas  for  the  head  and  shank  were  arbitrarily 
selected  while  the  other  sections  were  the  same  in  cross  section  as 
the  ceramic. 

The  weight  factors  were  assigned  as  given  in  Table  6.  This  is 
the  case  of  "equal"  weights,  meaning  that  the  terms  in  the  achievement 
function  are  all  approximately  of  the  same  order  of  magnitude.  The 
step  sizes  EPS  were  chosen  to  be  1%  of  the  initial  value.  Some 
consideration  was  given  to  the  acceleration  factor  ALPHA.  If  the 


-  •  -  - - 


«H _ mm 


» _ _ _ _ _ _ _ 
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Table  7 

Variable  Assignment  in  "VARIB"  for  the  Problem  Example 

SUBROUTINE  VARIB 
REAL  L(9) 

DIMENSION  RH0(9) ,  C(9),  A(9) ,  ETA(9) ,  X(12) 
CO!CfON/COMM  12/RHO,  C,  L,  A,  ETA 
COMMON /COMM  02/X 
C  User  Supplied  Variables  Are  Listed  Here 
L(6)  =  X(l) 

L(5)  =  X(2) 

L(3)  =  X(3) 

L(2)  =  X(4) 

L(l)  =  X(5) 

RETURN 

END 


Table 


ALPHA,  BETA  1.5 


(Continued) 
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acceleration  is  chosen  to  be  too  large,  it  causes  the  system  to  fluc¬ 
tuate  back  and  forth  about  an  optimal  solution.  In  this  condition, 
the  program  is  said  to  "hunt  to  and  fro"  and  the  results  are  likely 
to  be  erratic.  After  some  trial  and  error,  a  value  of  1.5  was 
selected  for  ALPHA.  The  achievement  tolerances  EPSY  were  both  set  to 
0.001. 

The  output  for  the  case  of  "equal"  weights  is  shown 
in  Table  9.  Each  iteration  is  presented  to  show  the  orderly 
improvement  in  the  objective  functions.  The  output  shows  both  the 
solution  vector  of  lengths  and  the  values  of  each  objective  function 
at  those  lengths.  An  examination  of  the  sequence  of  outputs  at  each 
step  reveals  the  resonance  frequency  climbing  towards  5000  Hz,  Q 

m 

being  reduced,  and  kgf^  being  increased,  weight  and  length  both  being 
reduced,  and  the  node  approaching  the  center.  To  achieve  this  state, 
some  of  the  variables  are  seen  to  vary  in  an  interesting  fashion  about 
the  center  variable.  The  tail  as  well  as  all  forward  sections  are 
reduced  in  length  while  the  fifth  section  insulator  is  increased. 

4.2  Variation  of  Weights 

The  next  set  of  tables  presents  the  outputs  after  the  weight  on 
one  or  another  deviation  variable  was  greatly  magnified.  This  tends 
to  steer  the  solution  to  favor  optimizing  that  function  over  any  other 
objective. 

Regarding  Table  10,  the  distorted  importance  of  causes  the 
objective  function  Y(3)  to  steadily  decrease.  To  achieve  this,  the 
length  of  the  head  shrinks  while  the  tail  grows.  Effectively,  the 


Table 


0.39170  0.39180  0.39150  0.39000 


Table  9  (Continued) 
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Table  10  (Continued) 
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water  load  is  brought  ever  closer  to  the  motor  section,  thereby 
increasing  its  effect.  Simultaneously,  the  mass  of  the  tail  is  in¬ 
creased  to  produce  a  larger  mass  loading  on  the  back,  thereby  effec¬ 
tively  clamping  the  back  face  of  the  ceramic.  These  conditions  com¬ 
bine  to  load  both  ends  of  the  transducer  and  to  lower  its  Q  .  Note 

m 

also  that  the  | admittance!  decreases  while  the  coupling  coefficient 
remains  nearly  stable. 

Next,  the  program  was  rerun  with  the  weight  factor  on  the  ad¬ 
mittance  deviation  increased  to  100,000,  in  an  attempt  to  produce  a 
larger  value.  (The  weight  on  was  decreased  to  its  former  value.) 
However,  because  of  the  water  loading  and  of  the  initially  small  size 
of  the  variable  perturbations  (step-size  EPS)  used  by  the  pattern 
search,  an  increase  in  the  magnitude  of  the  admittance  could  not 
be  effected.  The  program  could  find  no  improvement  over  the  initial 
calculated  value  of  0.15  my.  Of  course,  the  step  reductions  were  not 
successful  (only  a  step  increase  would  have  helped);  therefore,  the 
program  executed  its  maximum  allowed  number  of  step  reductions  and  sub¬ 
mitted  as  a  best  guess  the  original  vectors  Y(I)  and  X(I). 

When  the  weight  factor  associated  with  the  coupling  coefficient 
3 

^eff  was  i-ncrease<l  to  10  ,  the  iterations  shown  in  Table  11  resulted. 
Here,  the  final  solution  is  not  distorted  very  much  from  that  for 
"equal"  weights.  This  may  be  because  the  strong  emphasis  placed  on 
the  coupling  steers  the  solution  in  the  same  basic  direction  as  the 
combination  of  the  other  objectives  does.  Interestingly  enough, 
even  with  such  a  high  weighting  to  obtain  the  maximum  ^e££»  the  value 
at  the  end  of  the  4th  iteration  is  smaller  than  the  1st.  This  is 


front  ceramic 


Table  11  (Continued) 


probably  caused  by  the  influence  of  the  other  objective  functions 
and  the  limited  number  of  searches. 
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A  last  example  was  set  in  which  the  weight  on  the  length  devia¬ 
tion  was  increased  to  1,000.  As  in  the  previous  case  for  k^^,  the 
goal  of  minimum  length  must  be  compatible  with  the  other  objectives 
because  through  four  searches  the  solution  was  identically  equal  to 
that  of  "equal"  weights. 

4 . 3  The  Question  of  Uniqueness 

There  are  three  separate  phenomena,  each  of  which  are  multimodal, 
which  obviate  the  possibility  of  a  unique  solution.  Firstly,  the 
piezoceramic  motor  section  of  the  transducer  has  an  infinity  of 
resonances  in  any  direction  of  vibration.  Secondly,  the  addition  of 
variable  masses  for  the  head  and  tail  provides  an  infinity  of  combina¬ 
tions  to  produce  any  specific  resonance  frequency.  Addition  of  other 
variable  masses  (insulators,  shank,  and  acoustic  window)  serve  to 
increase  the  infinity  of  combinations  which  have  the  same  resonance 
frequency.  While  the  resonance  frequency  isn't  the  only  objective 
function,  its  deviation  variables  will  dominate  the  achievement  func¬ 
tion  in  the  first  few  iterations.  Thirdly,  any  nonlinear  optimization 
technique  is  multimodal  and  the  achievement  function  reached  in  the 
NLGP  technique  is  typically  a  function  of  the  starting  point,  the  step 
sizes,  and  the  order  of  ranking  of  the  independent  variables. 


CHAPTER  5 


SUMMARY  AND  CONCLUSIONS 


5.1  Summary 

The  broad  objective  of  this  thesis  was  to  develop  an  accurate 
tool  for  transducer  design.  A  modified  Mason's  Equivalent  Circuit 
was  used  to  construct  a  mathematical  model  (Figure  7)  which 
suitably  predicts  the  design  parameters  of  resonance  frequency, 
admittance  magnitude  and  phase,  and  node  placement.  Nonlinear  Goal 
Programming  was  selected  as  a  viable  multiobjective  optimization 
technique  to  evaluate  the  mathematical  model,  compare  the  results  with 
the  specified  design  parameters,  and  provide  the  optimal  compromise 
values  for  the  design  variables. 

In  Chapter  1,  the  basic  transducer  design  problem  was  presented. 
Transducer  theory  was  explored  in  Chapter  2  where  equations  which 
relate  the  design  variables  to  the  transducer  responses  were  developed. 
The  reliability  of  the  equations  and  of  the  computer  algorithm  was 
examined  with  comparisons  between  both  theory  and  experiment. 

Chapter  3  contained  a  description  of  Nonlinear  Goal  Programming  and 
its  application  to  the  transducer  model.  A  hypothetical  transducer 
was  designed  in  Chapter  4  to  illustrate  NLGP  operation  and  to  display 
the  importance  of  the  weight  factor. 


5.2  Conclusions 


Goal  programming  is  a  powerful  technique  for  obtaining  the 
optimal  compromise  solution  to  complex  multiobjective  problems. 
"Tonpilz"  transducer  design,  involving  all  segments  of  the  body  of  the 
transducer  and  its  associated  housing,  can  easily  be  accomplished 
using  this  technique,  but  it  must  be  recognized  that  as  with  any  other 
nonlinear  optimization  tool,  or  any  multimodal  system,  the  resultant 
solution  may  well  be  only  a  local  optimal  solution.  Further  research 
should  be  done  to  create  a  theoretical  model  to  explain  the  frequency 
differences  encountered  by  actual  "identical"  transducers,  which 
caused  the  errors  outlined  in  Chapter  3. 

One  of  the  more  exciting  aspects  of  this  design  program  is  that 
it  provides  an  opportunity  for  theoretical  research  into  the  behavior 
of  all  types  of  transducers.  The  optimization  technique  allows  any 
combination  of  attributes  to  be  examined  in  relation  to  the  problem 
variables  as  has  never  before  been  possible.  Finally,  this  technique 
can  easily  be  applied  to  the  design  and  analysis  of  other  types  of 
systems  and  to  other  optimization  problems.  It  is  hoped  that  this 
thesis  might  guide  the  researcher  in  other  areas. 


APPENDIX  A 


CALCULATION  OF  WAVE-GUIDE  EQUATION  (8)  FROM  LOADED  T-SECTION 


Figure  13  is  a  T  circuit  representation  of  a  wave  guide 
terminated  with  a  load  impedance  Z  .  The  input  impedance  can  be  seen 

Li 

to  be: 


(Z  +  jZ  tan(x/2))  (— jZ  /sin(x)) 

Zin  =  ZL  +  jZQtan(x/2)  -  jZQ/sin(x)  +  ^Z0tan^2^  • 


(20) 


where  x  =  k. i. . 

i  i 

Using  a  trigonometric  substitution. 


tan(x/2)  -  l/sin(x)  =  -cot(x) 


(21) 


'in 


-jZ^Z^cot (x)  +  Z^tan(x/2)  (2/sin(x)  -  tan(x/2)) 
ZL  -  jZQcot(x) 


(22) 


but 


2/sin(x)  -  tan(x/2)  =  l/tan(x/2)  . 


(23) 


Therefore , 

-jZ0ZLCOt(x)  +  Z0 

Zin  ZL  -  jZQcot(x)  ’ 


(24) 


which  is  expressed  in  the  text  as: 

ZQ(ZL  +  jZQtan(x)) 
Zin  *  Zq  +  jZ^tan(x) 


(25) 


APPENDIX  B 


SOLUTION  OF  TRANSDUCER  CIRCUIT  WHICH  INCLUDES  A  STRESS  BOLT 

The  method  of  analysis  used  to  simplify  the  transducer  circuit 
is  a  tt-T  transformation  as  follows.  Referring  to  Figure  14,  one  can 
equate: 

5^2  V3  hh 

^A  “  +  52  +  ^3  ?B  =  q  +  C2  +  £3  ^C  £1  +  £2  +  £3  (26) 

where  the  £'s  are  impedances. 

The  circuit  shown  in  Figure  7  can  be  succinctly  reduced  to 
Figure  15.  The  elements  dependent  on  1/sin (x)  are  designated  by  number 
as  assigned  in  Figure  12,  and  the  combinations  of  elements  dependent 
on  tan(x/2)  are  designated  by  letter.  The  lettered  circuit  elements 
of  Figure  15  are  defined  as  follows: 

ZA  =  Z[^head1  +  jz7tan(k7V2)  +  jz3tan(k3V2) 

ZB  =  jZ3tan(k3£3/2)  +  jZ4tan(k4H4/2) 

Zc  =  jZ4tan(k4£4/2)  +  jZ5tan(k5^5/2)  (27) 

ZQ  =  jZ5tan(k5£5/2)  +  jZ6tan(k6*6/2) 

ZE  =  jZ6tan(k6Jl6/2)  +  jZytanCk^/2) 


Redrawing  the  circuit  from  Figure  15,  one  obtains  Figure  16(a). 
Applying  the  ir-T  transformation  to  the  circled  segment  yields 
Figure  16(b).  The  newly  defined  elements  from  this  operation  are: 
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Reduced  Transducer  Circuit  as  a  Result  of  One  Application  of  a  ir-T  Transformation 
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The  circuit  has  now  been  reduced  to  that  of  Figure  18.  From  this,  one 
can  calculate  the  mechanical  impedance  which  is: 


Z  =  Z  +  -Zc2ZBl  ,  z 
input  Z4  +  Z„„  +  Z„,  +  ZH3 


C2 


"B1 


(31) 


All  circuit  elements  dependent  on  l/sin(x)  are  given  by  number,  and 

these  terms  do  not  cancel  throughout  Equation  (31) .  It  is  evident 

that  one  cannot  let  x  =  0  in  any  section  without  generating  a  division 

by  zero  in  Z, 

input 


Figure  18.  Final  Reduction  of  Circuit. 


APPENDIX  C 


PROGRAM 

C  NLCP  TRANSDUCER  DESIGN  *  JULY,  1979 

C  THIS  PROCRAM  IS  FOR  TRANSDUCER  DESIGN.  ITS  PURPOSE  IS  TO 

C  PROVIDE  A  METHOD  FOR  DESICNINC  A  SPECIFIC  TONPILZ  TRANSDUCER  CIVEN 
C  THE  DESIRED  SPECIFICATIONS.  THE  ANALYST  SUPPLIES  THE  RESPONSES 
C  DESIRED  OF  THE  TRANSDUCER,  THE  VARIABLES,  AND  AN  INITIAL  TRIAL 
C  SOLUTION.  THE  PROBLEM  IS  THEN  SOLVED  VIA  NONLINEAR  COAL  PROCRAMMING 
C  AS  DEVELOPED  BY  JAMES  F.  ICNIZIO  (SEE  TEXT-  COAL  PROGRAMMING  AND 
C  EXTENSIONS,  FOR  DETAILS). 

c 

c 

C  THE  OBJECTIVE  FUNCTIONS  FOUND  IN  THE  PROCRAM,  UITH  THE  SUCCESTED 


c 

WEIGHTS  FOR  EQUAL  EMPHASIS  ARE 

: 

c 

I RCW 

-I  Y ( I ) 

RHS(I)  SICN(I) 

WEIGHT ( t) 

c 

PRIORITY  LEVEL  l 

,MP  R I- 1 

c 

10 

LBACK/LCERAMIC 

1 

+  1 

I 

c 

c 

11 

LFRON I/LCERAHIC 

1 

+  1 

1 

c 

PRIORITY  LEVEL  2 

, KPR I -2 

c 

1 

admittance  resonance 

ENTER  DESIRED 

+  .-1 

0.1 

c 

VALUE  (HZ) 

c 

2 

LFROST/LBACK 

1 

+  .-1 

10 

c 

3 

QUALITY  FACTOR,  QM 

0 

+  1 

1 

c 

A 

COUPLING  COEFFICIENT 

1 

-l 

10 

c 

5 

TOTAL  LENCTH 

0 

+  1 

10 

c 

6 

TOTAL  WEIGHT 

0 

+1 

1 

c 

7 

ABS(ADMITTANCE)  AT  FY 

AS  DESIP.ED 

-1 

1000 

c 

VALUE  (MHOS) 

c 

E 

PHASE  ANCLE  AT  FY 

0 

♦l 

100 

c 

9 

TOTAL  COST 

20Z  OF  TRIAL 

+1 

1 

c 

c 

. 

SOLUTION 

c 

C 

VARIABLE  IDENTIFICATION: 

USER  IS  PROVIDED 

WITH 

A  CHOICE 

c 

OF  AT  MOST  12  VARIABLES,  X(I) 

FROM  THE  RHO.C.L, 

A, AND 

ETA  OF 

C 

EACH 

WAVE  GUIDE  SECTION.  ALL 

OTHER  PARAMETERS 

ARE 

FIXED  AT 

c 

INITIAL  DATA  POINTS.  USER  MUS 

T  PROVIDE  SPECIFIC  VARIABLE 

C 

ASSIGNMENTS  IN  SUBROUTINE  VARIB  AT  THE  END  OF  THE  DECK. 

c 

(EC. 

TO  OPTIMISE  THE  LENCTH  OF 

THE  TAIL,  ONE  MUST  INTRODUCE  A 

C  CARD  L ( 6 )  -  X ( 1 )  IN  THE  SUBROUTINE  VAR I B . 

C 

C  WAVE  GU IDE  SECTION  DESCRIPTION 
C  1  HEAD 

C  2  SHANK  OF  HEAD 

C  3  INSULATOR 

C  A  CERAMIC  STACK 

C  5  INSULATOR 

C  6  TAIL 

C  7  RUBBER  WINDOW 

C  8  PRESSURE  RELIEF 

c 

C  NOTE:  TO  ELIMINATE  ANY  SECTION,  SET  THE  LENCTH  -0.0  IN  THE 
C  DATA  INPUT.  THE  CORRESPONDING  RHO,C,  AND  A  CAN  EQUAL  1.0 

c 


“IS  PAGi  Is  B*3I  QUALITY  PIUCIIUJM4 
raOM  OOP  Y  PWKVAiollKu  BO  DDC  _ _  . 


c* 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

V 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c* 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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VARIABLES 
NVAR 
NOB  J 
NPRIOR 
NACH 
NHAX 
X 

EPS 

RHS 

HP R  I 

UEICHT 

ALPHA 

BETA 

MAXRED 

IPRINT 

EPS  Y 
S1CN 


IROW 

Y 

XBASE 

ZTEST 

RHO 

C 

L 

A 

ETA 

TURNS 

CAP 

ZU 

ZPR 

STEP 

COST 


UTILITY 

DECISION  VARIABLES 
OBJECTIVE  FUNCTIONS 
PRIORITY  LEVELS 
TERMS  IN  ACHIEVEMENT  FUNCTION 


‘•••••••••••USERS 

VARIABLE 
CARD  1 

NVAR, NOB J .NPRIOR, NACH, NMAX 
CARD  2 

IPRINT, MAXRED 
CARD(S)  3- 

X(I)  .1-1, NVAR 
CARDS ( S )  4- 

EPS(I) ,1-1 , NVAR 
CARD(S)  5- 

RHS(I) ,1-l.NOBJ 
CARD  6 

ALPHA, BETA 
CARD  1 


NUMBER  OF 
NUMBER  OF 
NUMBER  OF 
NUMBER  OF 

MAXIHUM  NUMBER  OF  PATTERN  SEARCH  CYCLES  ALLOWED 
VECTOR  OF  DECISION  VARIABLES  AT  TEST  POINT 
VECTOR  OF  STEP  SIZES  USED  IN  PATTERN  SEARCH 
RICHT  HAND  SIDE  VALUES  OF  OBJECTIVE  FUNCTIONS 
PRIORITY  LEVEL  OF  GIVEN  DEVIATION  VARIABLE 
WEIGHT  FACTOR  OF  CIVEN  DEVIATION  VARIABLE 
PATTERN  SEARCH  ACCELERATION  FACTOR 
STEP  REDUCTION  FACTOR 

MAXIMUM  NUMBER  OF  STEP  REDUCTIONS  ALLOWED 
OUTPUT  SWITCH — 0  GIVES  ONLY  FINAL  RESULT 
1  GIVES  RESULT  OF  EVERY  PATTERN  SEARCH  CYCLE 
VECTOR  OF  IMPROVEMENT  TOLERANCES  IN  ACHIEVEMENT 
INDICATES  POSITIVE  OR  NEGATIVE  DEVIATION 
TO  BE  MINIMIZED  FROM  OBJECTIVE  FUNCTION  IROW  AT 
MPRI  LEVEL  WITH  WEIGHT  FACTOR  WEIGHT 
OBJECTIVE  FUNCTION  OF  GIVEN  DEVIATION  VARIABLE 
VECTOR  OF  FUNCTIONS  OF  DECISION  VARIABLES 
BASE  POINT  FOR  HOOKE-JEEVES  ALGORITHM 
ACHIEVEMENT  FUNCTION  VALUE  VECTOR  AT  A  TEST  POINT 
VECTOR  OF  DENSITY  VALUES  (MSS)  KG/CU  M 
VECTOR  OF  SOUND  SPEEDS  M/SEC 
LENGTH  OK  EACH  WAVE  GUIDE  SECTION,  M 
AREA  Or  EACH  SECTION  SQUARE  METERS 
VECTOR  OF  LOSS  FACTORS 

THE  SQUARE  OF  THE  TRANSFORMATION  FACTOR 

THE  CALCULATED  CAPACITANCE 

THE  LOADING  PRESENTED  TO  THE  FRONT  FACE 

THE  LOADING  PRESENTED  BY  THE  WEB  APPLIED  TO  THE 

PRESSURE  RELIEF 

.  THE  FREQUENCY  INCREMENTS  USED  TO  FIND  RESONANCE 
BY  ITERATION 

VECTOR  OF  VALUES  OF  EACH  COST  PER  UNIT  VOULME 
CUIDE  TO  INPUT  D AT A *»•••»••**••••***•******»•*•**•• • 

FORMAT 


513 

213 


8E10.4 


BE  10 • 4 
8E10.4 


2E10.4 


BUS  PAui  Lo 

IR0J1  (XtfX  DO  iiPfl  - 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


1 


2 

600 

601 

602 


EPSY(I) ,1-1 .NPRIOR 
CARDS  8-22 

SIGN(t) , I  ROW ( I ) ,MPRI(I) ,UEICHT(I) ,I-t ,NACH 
CAROS  23-24 
COST ( I )  .1-1,8 
CARDS  25-34 

RHO(I)  ,C(I)  ,1(1) ,A(I) ,ETA(t)  ,1-1,8 
CARD  35 

TURNS, CAP, ZW.ZPR 

NOTE:  ZW  AND  ZPR  ARE  COMPLEX  AND  REQUIRE  2 
CARD  36 
STEP 


8E10.6 

PS. 0,215, F10. 0 
8E10.6 
8E10.6 
6E10.6 

DATA  FIELDS  EACU 
15 


INFORMATION  ON  DIMENSION  STATEMENTS 

X(I).  I  -  TOTAL  NUMBER  OF  VARIABLES  (NVAP.)  . 

Z(I),  I  -  TOTAL  NUMBER  OF  PRIORITIES  (NPRIOR) 

EPS ( I )  ,  I  -  TOTAL  NUMBER  OF  VARIABLES  ( N VAR) . 

RHS(I),  I  -  TOTAL  NUMBER  OF  OBJECTIVES  (NOBJ) 

EPST(I),  I  -  TOTAL  NUMBER  OF  PRIORITIES  (NPRIOR) 

SICN(I),  I  -  TOTAL  NUMBER  OF  TERMS  IN  ACHIEVEMENT  FUNCTION  (NACH) 

IROW(I),  I  -  TOTAL  NUMBER  OF  TERMS  IN  ACHIEVEMENT  FUNCTION  (NACH). 

KPRI(I),  I  -  TOTAL  NUMBER  OF  TERMS  IN  ACHIEVEMENT  FUNCTION  (NACH) 

WEIGHT ( I ) ,  I  -  TOTAL  NUMBER  OF  TERMS  IN  ACHIEVEMENT 
FUNCTION  (NACH). 

Y ( I ) ,  I  -  TOTAL  NUMBER  OF  OBJECTIVES  (NOBJ). 


DIMENSION  X ( 1 2 )  ,Y(12) 

COMMON /COMMOl/N VAR, NOBJ, NPRIOR, NACH, NMAX 
COMMON /COMMO 2 /X 
COMMON /COMMO  7 / Y 
CALL  DATAIN 

ALL  NECESSARY  DATA  IS  INPUT  AND  PRINTED  OUT 
CALL  HJALC 

THE  PATTERN  SEARCH  ALGORITHM 
PRINT  600 
PRINT  601 
CALL  VARIB 

CALL  Y VALUE  * 

DO  1  K-l.NOBJ 
PRINT  602 ,K, Y (K) 

CONTINUE 
PRINT  603 
DO  2  X- 1  , N VAR 
PRINT  604 , K , X(K) 

CONTINUE 

FORMAT! IHO ,25X, 'MAXIMUM  NUMBER  OF  PATTERN  SEARCH  CYCLES  EXCEEDED', 
J/.1H  ,25X,'BEST  SOLUTION  TO  THIS  POINT  FOLLOWS') 

F0KMA7( 1H0.25X,  '  OBJECTIVE  FUNCTION  VALUES') 

FORMATUH  ,25X,'Y(*  ,12,')  -  '  ,  F  l  3 . 6) 


IHIS  PAGH  IS  BiSI  QUALITY  PfUCTKMBIil 

KlttV  COr  Y  TO  DDC  _ 


_ - - -  . . 


77 


* 


603  FORMAT(1HO,25X, 'DECISION  VARIABLES  VALUES') 

604  FORMAT ( 1 H  ,25X,'X(',I2,')  -  ',F13.6) 

STOP  • 

END 

SUBROUTINE  OATAIH 
COMPLEX  ZW.ZPR 
REAL  L ( 9 ) 

INTECER  STEP 

DIMENSION  X( l 2) , EPS (12)  ,  RHS ( 1 2 )  .  EPS Y ( 3 )  ,  S ICN ( 2 4 ) , I  ROM 
$(24) ,MPRI(24) .WEIGHT (24)  ,RH0(9)  ,  C ( 9 )  ,  A ( 9 ) . ETA ( 9 ) .COST (9) 
COMMON /C0MM01/N VAR, NO BJ , N PR IOR , N ACH . NMAX 
COMMON / C0MM02 / X 
COMMON /COMMO 3 /EPS 
COMMON /COMMO 4 /RHS 

COMMON /COMMO 5 /SIGN, ME ICHT , IROW.MPRI 

COMMON /COMMO 5 /EPSY 

COMMON /COMM  1 2 /RHO.C.L, A, ETA 

COMMO N /COMM I  3 /TURNS .CAP, ZM.ZPR, STEP, FT, COST 

COMMON /COMMO 9 /ALPHA, BETA 

COMMON /COMM 10/MAXRED 

COMMON /COMM 1 1 /I  PRINT 

READ  500, N VAR, NOB J . NPR IOR, NACH, NMAX 

PRINT  600,  NVAR.NOBJ .SPRIOR, NACH, NMAX 

READ  500.IPRINT.MAXRED 

PRINT  612,IPRINT, MAXRED 

READ  501, (X(I)  , I- l , N  VAR) 

PRINT  602 
DO  2  I - 1 , N VAR 
PRINT  60I,I,X(I) 

2  CONTINUE  ' 

READ  50I,(EPS(I),I-l,NVAR) 

PRINT  603 

DO  3  I-t.NVAR 
PRINT  604,I,EPS(I) 

3  CONTINUE 
PRINT  800 

READ  50l,(RHS(I) ,I*1,N0BJ) 

PRINT  605 
DO  4  1*1, NOB J 
PRINT  606,1,RHS(I) 

4  CONTINUE 
FY-  RHS(l) 

READ  501,  ALPHA, BETA 
PRINT  607, ALPHA, BETA 
READ  501,(EPSY(I) ,I-1,NPRI0R) 

PRINT  608 

DO  5  1*1, NPR IOR 

PRINT  609,I,EPSY(I) 

5  CONTINUE 
PRINT  610 

DO  1  N- l , NACH 


I  I 

1  J 
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READ  502,SICN!N) ,IROW(N)  ,MPRI(N)  . WEIGHT (H) 

PRINT  6  1  1  ,SICN(N)  , IROU(N)  ,  MPR  I  {  N  )  .WEICIIT(N) 
t  CONTINUE 
PRINT  701 

READ  501,  (COST(I) ,1-1 ,8) 

DO  20  N-1,8 

READ  50l,RHO(N) , C(N ) ,  L  (N)  ,  A(N) , ETA(N) 

PRINT  700,N,RHO(K) , C (N ) , L (N ) , A ( N ) , ETA ( N ) ,COST(N) 

20  CONTINUE 
PRINT  702 

READ  501, TURNS, CAP, ZW.ZPR 
PRINT  703, TURNS, CAP, ZW.ZPR 
READ  500, STEP 
PRIST  704 ,  STEP 

500  FORMAT (1615) 

501  FORMAT ( 8E 1 0*  4) 

502  FORMAT (F5.0,215,F10.0) 

700  FORMAT(15X,I2,6E15.4) 

70  1  FORMAT!  111!  ,  22X,  'RHO'  ,  15X,  'C' , 1 4 X  ,  * L  '  , 1 4 X , * A  * , 13 X , * ET A' , 13X, 'COST') 

702  FORMAT!/, 25X,'  TURNS  R AT l 0 ' , 9 X , ' C AP ' , 7 X , ' RAD  LO AD ' , 1 5X , ' WE B  LOAD') 

703  FORMAT(25X.E10.4,5X,E10.4,4X,2E10.4,3X,2E10.4) 

704  FORMAT!/, 25X,'  FREQUENCY  LOOP  STEP  S  IZ E  '  , 1 1 0 , / //) 

600  FORMAT ( l H  ,25X, 'NUMBER  OF  DECISION  VARIABLES  - '  ,  T6 8 , 1  5 , / , 1 H  ,2SX,' 
SSUMBER  OF  OBJECTIVE  FUNCTIONS  - ' , T 6 8 , 1  5 , / , 1 H  ,25X, 'NUMBER  OF  PRIOR 
SITY  LEVELS  - ' , T6 8 , 1  5 , / , 1 H  ,25X, 'NUMBER  OF  TERMS  IN  ACHIEVEMENT  FUN 
SCTION  -' ,T68,15, /, 1H  , 25 X ,' MAXIMUM  NUMBER  OF  SEARCH  PATTERNS  -',T6 
$3.15,/) 

612  FORMAT ( 1 H  ,'  IPRINT '  ,  I  5 , 25X, '  MAXIMUM  NUMBER  REDUCT ION S ' , T63 , 1 5) 

601  FORMAT ! 1 H  ,  2 5X  ,  ' X ! '  ,  I  2  ,  '  )  -  '.E10.4) 

602  FORMAT!lHl,25X, 'VECTOR  OF  INITIAL  CUESSES  OF  DECISION  VARABLES  — 
$X') 

603  FORMATOH1.25X, 'INITIAL  INCREMENTS  FOR  DECISION  VARIABLES  —  EPS') 

604  FORMAT ! 1 H  , 25X , ' EPS (  '  .  12  ,  '  )  -  ',F10.6) 

605  FORMAT! 1H1 ,25X, 'VALUE  OF  RIGHT-HAND  SIDE  OF  OBJECTIVE  FUNCTIONS') 

606  FORMAT  (  1  It  ,  2  SX ,  '  RHS  !  '  ,  I  2  ,  '  )  -  '.F10.4) 

607  FORMAT!  1111 ,25X,  'ALPHA  -  '  ,  F  8 . 4 , 1  OX  , '  B  ET  A  -  ',F8.4) 

608  F0RMATIIH1 ,25X. 'VECTOR  OF  ACHIEVEMENT  FUNCTION  TOLERANCES') 

609  FORMAT ! 1 H  ,  25X  , ' EPS Y  !  '  ,  I  2  .  '  )  -  '.F9.7) 

610  FORMAT! 1111 ,T25, 'DEVIATION  T YP E '  ,T4 5 , ' 0 B J ECTI VE  NUM BER ' ,T65 ,' PRIORI 
STY  LEVEL' ,TB5 , 'WEIGHT  FACTOR') 

611  FORMAT ! 1 H  ,T30. F6. 1 ,T54, 12 ,T74, I2.T91 ,F8,2) 

8C0  FORMAT! IHl , 25X. 'OBJECTIVE  FUNCT  IONS '  ,  /  ,  2 5 X , ' Y { 1 )  -  F Y ' , / ,  2 5 X , ' Y ( 2 ) 

1  -  LBACK/LFRONT' , / , 2SX, ' Y( 3)  -  Q Y '  ,  /  ,  2 5X , ' Y < 4 )  -  KEFF '  ,  /  ,  2 5 X , ' Y ( 5 ) 

2  -  LENCTH' , / , 2  5  X , 'Y!6)  -  WEI CH T ' , / , 2 5 X , ' Y ! 7 )  -  ABS ! ADM  I T) ' . / , 2 5X , 

3  '  Y  <  8 )  -  PHASE  OF  ADM  I T  '  ,  /  ,  2 5 X ,  '  Y ( 9 )  -  COST ' , / , 2 5X , ' Y ( 1 0 )  -  LBACK/ 
4LCERAMIC' ,/,25X, ' Y( l  1)  -  L FRONT /LC ERAMIC ' ) 

RETURN 

END 

SUBROUTINE  VALUE(Z) 

C  SUBROUTINE  VALUE  EVALUATES  THE  ACHIEVEMENT  FUNCTION  CIVEN  THE 

C  PREVIOUSLY  EVALUATED  FUNCTIONS  OF  THE  DECISION  VARIABLES 

DIMENSION  Z( 3)  ,Y(  12)  ,RII5(  1  2)  ,  S  I CN ( 2 4 )  ,  l ROW ( 2 4 ) , MPR I ! 2  4 ) ,WEICHT(24) 
COMMON  /COMHO  l  /N  VAR,  NOBJ  ,  N  PR  IOR  ,  N  ACII ,  NMAX 


fHI3  PAO£  IS  I  ^AUll 

KiOSI  (XAr  i  Pifcvrfis ti£D  DO  DDC  _ 
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5 
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C 

C 

c 

c 

c 

c 

100 


30 

101 


1 

c 

199 
C 

c 

200 


C0MM0N/C0MM04/RHS 

COMMON /COMMO 5 / S ICN . WEI CHT , I ROW ,MPRI 

COMMON/ COMMO  7 /Y 

CALL  VAR  I  B 

CALL  Y VALUE 

DO  1  K-l.NPRIOR 

Z ( K) -0 . 0 

CONTINUE 

DO  2  K-l.NACH 

IF(S1CN(K) . LT . 0 . 0 ) C0T04 

DE  V-Y (1R0U(K))-RHS(IR0V(K)) 

C0T05 

DEV-RHS(IROW(K) ) -Y(IROW(K)) 

IF(DEV.LT.0.0)DEV-0.0 

Z(MPRI(K) )-Z(MPRI(K) )+WElGHT(K) *DEV 

CONTINUE 

RETURN 

END 

SUBROUTINE  HJALC 
DIMENSION  XSAVE(12)  ,Y(12) 

DIMENSION  X(12) ,XBASE(12) 

DIMENSION  EPS (12) ,Z( 3) 

COMMON /COMMO 1 /N VAR, NOB J , NPR 10 R , N ACH , NMAX 

COMMON /COMMO 2 /X 

COMMON /COMMO 3 /EPS 

COMMON / COMMO  2 / Y 

COMMON /COMMO 8 / Z 

COMMON /C0MM09 /ALPHA, BETA 

COMMON /COMM  10 /MAXR ED 

COMMON /COM Ml  1 /IPRINT  1  * 

INITIALIZATION  OF  PATTERN  SEARCH 
HCOUNT-O 

KCOUNT  COUNTS  THE  NUMBER  OF  PATTERN  SEARCHES  ATTEMPTED 
NREDUC-0 

NREDUC  COUNTS  THE  NUMBER  OF  STEP  SIZE  REDUCTIONS 
KFLAC-1 

KFLAC  IS  A  SWITCH  WHICH  TELLS  THE  PROCRAM  WHETHER  THE  NEW  TEST 
POINT  FOUND  BY  ACCELERATION  IMPROVES  THE  VALUE  OF  THE  ACHIEVE¬ 
MENT  FUNCTION 
CALL  VALUE(Z) 

PRINT  601 
DO  30  K-l.NOBJ 
PRINT  602,X,Y(K) 

CONTINUE 
PRINT  603 
DO  l  K-I.NVAR 
XBASE (K) -X (K) 

PRINT  i04,K,X(X) 

CONTINUE 

START  THE  SEARCH  PATTERN 
NIJK-0 

NIJK  COUNTS  THE  NUMBER  OF  IMPRO VEM ENTS , I T  ANY.  DURINC  A  PATTERN 
SEARCH 

DO  3  K-l.KVAR 


rrNRBiH— i  — - -  » 
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C  STEP  FORWARD  AND  EVALUATE 

X(K)-X(K)+EPS(K) 

C  SUBROUTINE  DECIOE  TAKES  THE  VALUE  OF  THE  TEST  POINT.  EVALUATES 

C  THE  AC  HI  EVF.MF.  NT  FUNCTION  VIA  ZVALUE.  AND  DECIDES  IF  THE  TEST 

C  POINT  IS  AN  IMPROVEMENT.  IJK-1  FOR  IMPROVEMENT,  IJK-0  IF  NOT 

CALL  DECIDE(IJK) 

NIJK-NIJK+IJK 
IF(IJK.EQ. 1)  CO  TO  3 
C  STEP  BACKWARDS  AND  EVALUATE 

X(K)-X(K)-2.0*EPS(K) 

CALL  DECIDE(IJK) 

NIJK-NIJK+IJK 
IF(IJK.EQ.l)  CO  TO  3 
X(K)-X(K)+EPS(K) 

3  CONTINUE 

MCOUNT-MCOUNT+1 

C  CHECK  IF  MAXIMUM  NUMBER  OF  SEARCH  CYCLES  EXHAUSTED 

I F ( MCOUNT  . EQ.NMAX)G0T0999 

C  CHECK  TO  SEE  IF  PATTERN  RESULTED  IN  NO  IMPROVEMENT 

IF(NIJK.EQ.O.AND.KFLAC.EQ. DCOT0300 
C  ACCELERATION 

DO  10  K-l.HVAR 
XSAVE (K) “X (K) 

10  CONTINUE 

DO  11  K-l.NVAR 

X(K) -X(K)+ALPHA*(X(K) -XBASE (K>) 

11  CONTINUE 

DO  12  K-l.NVAR 
XBASE (K) -XSAVE(K) 

12  CONTINUE 

C  NOW  CHECK  VALUE  OF  ACHIEVEMENT  FUNCTION  AT  ACCELERATION  POINT 

KFLAC-0 

CALL  DECIDE! I JK) 

1F(1JK.EQ.0)KFLAC-1 

IF ( I  PR  I NT . NE . 0) PRINT  600, MCOUNT 

IF( KFLAC . EQ .  1)  CO  TO  1 7 

PRINT  601 

DO  1 4  K-l.NOBJ 

PRINT  602 ,K , Y (K) 

1*  CONTINUE 

PRINT  603 
DO  IS  K-l.NVAR 
PRINT  604  ,K,X(X) 

IS  CONTINUE 

PRINT  606 
CO  TO  199 
17  CONTINUE 

PRINT  607 
13  C0T0199 

C  REDUCE  STEP  SIZB 

300  NREDUC-NREDUC+1 

C  IF  MAXIMUM  NUMBER  OF  STEP  SIZE  REDUCTIONS  EQUALED,  TERMINATE 

DO  301  K-l.NVAR 
X(K) -XBASE(K) 


ISIS  PAUI  is  »A3T  QUALITY  FRACT10ABL1 
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301  CONTINUE 

IF(NREDUC.CT.KAXRED)  CO  TO  316 

304  DO  305  K*  1  ,  N  VAR 
EPS(K)-EPS(IC)  »BETA 

305  CONTINUE 

PRINT  605,NREDUC,MCOUNT 
C0T0199 

316  PRINT  608 
PRINT  801 
CALL  VARIB 
CALL  Y VALUE 

DO  318  K-l.NOBJ 
PRINT  602,K,Y!K) 

318  CONTINUE 

PRINT  603  y 

DO  319  K-l.NVAR 
PRINT  604,K,XBASE(K) 

319  CONTINUE 

317  STOP 

600  FORMAT! IH1 ,25X,  13  ,  *  SEARCH  PATTERNS  HAVE  BEEN  COMPLETED') 

601  FORMATdlll  ,25X,#  OBJECTIVE  FUNCTION') 

602  FORMAT  (111  ,  2  5X  .  '  Y  (  '  ,  I  2  ,  ' )  -  '.F13.6) 

603  F0RXAT(1H1,25X, 'DECISION  VARIABLES') 

604  FORMAT ( 1 H  , 25X , ' X ( '  ,  1 2  ,  ' )  -  '.F13.6) 

605  FORMAT(  nil , 25X, ' STEP  REDUCTION  NUMBER  ',13,'  OCCURRED  AFTER  SEARCH 
$  PATTERN  NUMBER  ',13) 

606  FORMAT ( 1 H  ,25X,'KFLAG  -  0  SO  THAT  THE  TEMPORARY  BASE  POINT  IMPROVE 
$S  THE  SOLUTION') 

607  FORMAT ( 1 H  ,25X,'KFLAG  -  1  SO  THAT  THE  TEMPORARY  BASE  POINT  DOES  NO 
•  ST  IMPROVE  THE  S OLUT 10 N ' , / , 1 H  ,25X, 'SHOULD  THE  PATTERN  SEARCH  FIND 

SNO  IMPROVEMENT,  STEP  SIZES  WILL  BE  REDUCED') 

608  FORMAT( 1H1 , 25X, 'PATTERN  SEARCH  TERMINATED  —  MAXIMUM  NUMBER  OF 
SSTEP  REDUCTIONS  EQUALED') 

999  RETURN 
END 


SUBROUTINE  DF.CIDE(M) 

C  SUBROUTINE  DECIDE  CALLS  SUBROUTINE  VALUE  TO  EVALUATE  THE  ACK- 

C  IEVEHF.NT  FUNCTION  AT  THE  TEST  POINT  AND  COMPARES  THIS  TO  THE 

C  PREVIOUS  VALUE.  IF  THE  POINT  WILL  PROVIDE  IMPROVEMENT,  M-l 

C  IS  RETURNED  TO  THE  ARCUMENT  LIST,  IF  NOT  M-0  IS  RETURNED. 

DIMENSION  X ! I  2) 

DIMENSION  Z ! 3 ) ,ZTEST(3) ,EPSY!3),Y!12) 

COMMON /COMMO 1 /N VAR, NOB J , N PR IOR , N ACH , NMAX 

COMMON / CO MM 02 /X 

COMMON /C0MM06/EPSY 

COMMON /COMMO  7 /Y 

COMMON /C0MM09/Z 

CALL  VALUE ( ZTEST) 

DO  10  KM  ,  N  PR  IOR 
IF(ZTESTO)  .GT.Z(1))G0  TO  12 
IF!ZTEST!K) .LT.Z(K))  CO  TO  11 
10  CONTINUE 
C0TO12 

11  DO  16  K-l.NPRIOR 

C  TERMINATION  BY  MEANS  OF  ZERO  FOR  ALL  PRIORITY  LEVELS 

IF( ZTEST |K)  .NE. 0.0) GOTO  17 
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16  CONTINUE 
C0T0999 

17  DO  13  K-l.NPRIOR 

C  TERMINATION  TEST  BY  TOLERANCES 

IF(ABS(Z(K)-ZTEST(K) ) .GT . EPS Y (K ) ) COTO l 4 

13  CONTINUE 
C0T0999 

14  DO  IS  K-l.NPRIOR 
Z ( K)  -ZT ES'C  ( K) 

15  CONTINUE 
M-l 

RETURN 
12  H-0 

RETURN 

999  PRINT  600 
PRINT  601 
00  993  K-l.NOBJ 
PRINT  602  ,K, Y (K) 

998  CONTINUE 

PRINT  603  * 

DO  997  K- 1 , N VAR 
PRINT  604 , X , X(K) 

997  CONTINUE 

STOP 

600  F0RMAT(1H1,2SX. 'SEARCH  ALGORITHM  IS  COMPLETED  —  TOLERANCES  FOR  IN 
$P  ROVED  SOLUTIONS  SATISFIED') 

601  FCRMAT( [HI , 25X, 'OPTIMAL  OBJECTIVE  FUNCTION  VALUES') 

602  FORMAT ( 1 H  . 25X  ,  ' Y (  '  ,  I  2 , ' )  -  '.F13.6) 

603  F0RMAT(1III,25X, 'OPTIMAL  DECISION  VARIABLE  VALUES') 

604  FORMAT ( 1H  , 25X . 'X (  '  ,  1 2 .  '  )  -  *,F13.6) 

END 

SUBROUTINE  YVALUE 
INTECER  FSTART, STEP, FREQ 

COMPLEX  ZW,ZPR,KAY(9) , ZLI NE , ZWEB , ZRAD , ZHE AD , ZMECH , CCOS , CSIN , ZLO AD . 
$J,S(3) ,LB,LF, ADMIT, ARC,KL(9) , ZT A IL , ZPR ESS , ZW AT ER , XTAN , XS IN 
REAL  L<9) ,LFRO NT, LBACK.SINH, COSH, TANK 

DIMENSION  RHO ( 9 ) ,C(9) ,A(9> ,ETA<9) ,F(3) ,Y(12) ,Z(9) .COST (9) 

COMMON /C0MM12/RH0.C.L, A, ETA 

COMMON /COMM  13 /TURN S, C AP ,ZW, ZPR, STEP, FY, COST 

COMMON  / COM’.tO  7  /Y 

SlMll(X)-  (EXP(X)-EXP(-X)) /2 

COSH ( X ) -  <EXP(X)+EXP(-X)) /2 

TAN  II  (X)-  (EXP(X)-EXP(-X)) /(EXP(X)+EXP(-X)) 

XTAN ( ARC) - (TAN ( REAL ( ARC) ) +J  *TANU ( AIMAG ( ARC) ) ) / ( 1 -J  *TAN (REAL( ARC) ) 
$*TANH( AIMAC( ARC) ) ) 

XS IN ( ARC) -S IN ( REAL( ARC) ) *COSH ( AIM AC ( ARC) ) +J * COS (R E AL ( ARC ) ) « S I NH ( 
SAIMAG(ARC)) 

ZLINE(ZL0AD , ZMAT , ARC) -  ZM AT* ( ZLO AD*J * ZMAT*XTAH ( ARC ) ) / ( ZMAT  +  J *ZLO 
$AD*XTAN(ARC) ) 

J-(0. 0,1.0) 

PI-  3.14159 
DO  3  1-1,8 

Z<I)-RHO(I)*C(I)*A(I) 

3  KAY(I)-(2*PI*L(I)*(1-J*ETA(I)))/C(t) 

ZWATER-ZW*A(l) 

ZPRESS-ZPR*A(8) 

FSTART-IFIX(FY)-500 
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S(l)-(0. 0,0.0) 

S(2)-(0.0,S.0E  04) 

S( J)-(0.0.-5.0E  04) 
r(i)-Fm-Fo>-o.o 

C  BECIN  FREQUENCY  LOOP 
18  FREQ -FSTART -STEF 

17  CONTINUE 

FREQ-FREQ+STEP 
DO  4  1-1,8 

4  KL(I)-FR£Q»KAY(I) 

ZTAIL-ZLINF.(ZLINE((0.0,0.0)  ,Z<6),KL(6))  ,Z(5)  ,KL(5)) 

ZWEB-  ZLINE(ZPRESS ,Z(3)  ,KL(8) ) 

ZRAD-ZLINE(ZLINE(ZUATER,Z( 7)  .KL(7)),Z(1),KL(1)) 

ZWEB-ZUEB+ZRAD 

ZHEAD-ZLINE(ZLINE(ZWEB,Z{2)  , XL( 2) )  ,  Z  (  3)  ,  KL(  3)  ) 

ZMECH-  -J*Z(4)/XSIN(KL(4))  + ( Z H EAD+J *Z ( 4 ) *  XTAN (KL ( 4 ) /2 ) ) * ( ZTA IL+4 


1*Z(4) *XTAN(KL(4) /2) ) /  (  ZHEAD+ZTA  IL+2>»J  *Z (4 )  *XTAJI (X L  t  * ; /*/ / 
ADMIT-  J*2*PI*FREQ»CAP  +TURNS /ZMECH 
C  COMPARE  VALUES  FOR  MAX, MIN  ON  REAL  AND  IMAGINARY  IMPEDANCE  AXES 
IFUEAL(ADMIT)  .GE.R£AL(S(1)))C0  TO  12 
IF(FREQ.EQ.FSTART+STEP)CO  TO  15 
CO  TO  13 

12  S ( 1 ) -ADMIT 
F( 1 ) -FREQ 
LF-ZHSAD 
LB-ZTAIL 

IF  (  AIM  AG  <  ADMIT)  .LE.AIMAC(SO)))  CO  TO  14 
S(3)-ADMIT 
F( 3) -FREQ 
CO  TO  17 

14  IF(FREQ.EQ.FSTART+STEP)  CO  TO  IS 
CO  TO  17 

13  I F ( A IMAC (ADMIT) .CE,AIMAC(S(2)))  CO  TO  19 
S ( 2 ) -ADMIT 

F( 2) -FREQ 
CO  TO  17 

15  FST  ART- FST  ART-50 
S(3)-(0.0,-5.0E04) 

S(l)-(0. 0,0.0) 

F ( 1 ) -F ( 3) -0  •  0 
CO  TO  18 
19  CONTINUE 

C  END  OF  FREQ  LOOP  NOW  VALUES  FOR  Y  ARE  CALCULATED 
CC-AIMAC(S(1) ) /(2*PI*F(l)) 

DY-CABS(S(2)-S(3>) 

IF(REAL(LB) . EQ .0, . AND. AIMAG(LB) . EQ .0.)  CO  TO  40 
LB ACK - ( L ( 4  ) /REAL(KL(4))) * AT AN ( AB S ( R EAL ( J *Z ( 4 ) /LB) ) ) 

CO  TO  41 

40  LB ACK-L { 4) /2 . 0 

41  IF(REAL(LF) . EQ .0. . AND . AIMAC (LF) .EQ.O.)  CO  TO  42 
LFR0NT-(L(4) /R EAL ( KL ( 4 ) ) ) * AT AN ( ABS ( REAL < J *Z ( 4 ) /LF))) 

CO  TO  43 

42  LFR0NT-L(4) /2.0 

43  CONTINUE 
Y(l)-F(l) 

Y(7)-LBACK/LFR0NT 

t(3)-F(l)/(F(2)-F{3)) 


c?  V 
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CK0T-0Y/(Y(3)  •2*PI«F(in 
Y  (4J-SQRT (ABS(CMOT/ C  CMOT+CC) ) ) 

Y(5)-L(l)+L(2)  +  L(3)  +  L(4)+L(5)+,L(6) 

T(6)-Z(  t ) *L( 1) /C ( 1)+Z(2) *1.(2) /C(2)+Z(3) *L<3) /C(l)+Z(4) »L(4) /C(4>+ 
IZ(5)*L(5)/C(5)+Z(6)*L(6)/C<6) 

Y(7)-CABS(S(l>  ) 

Y(8)-ATAN(AIMAG(S(1))  ./REAL(S(l))) 

Y( 10)-LBACK/L(4) 

Y(ll)-LFRO!IT/L(4)  - - 

Y(9)-0.0 
DO  20  1-1,8 

20  Y(9)-Y<9)  +  COST(I)*A(I)*L(I) 

TY-  F(3)  +  500. 

RETURN 

END 


SUBROUTINE  VARIB 
REAL  L ( 9 ) 

DIMENSION  RHO ( 9) ,C(9) ,  A(  9  ) ,ETA(9) ,2(12) 
COMMON /COMM l 2 / RHO , C , L , A, ET A 
COMMON /C0HM02/X 

C  USER  SUPPLIED  VARIABLES  ARE  LISTED  HERE 
L(6)«X(l) 

L(5)-X(2) 

L(3)-X(3) 

L(2)-X(4> 

.  L(l)-X(5) 

RETURN 

END 


BIBLIOGRAPHY 


1.  Ignizio,  J.  P.  Goal  Programming  and  Extensions.  D.  C.  Heath 

&  Co.,  Lexington,  MA,  Ch.  6  (1976). 

2.  Piezoelectric  Ceramics — Ferroxcube  Corp.  Application  Book. 

J.  Van  Randeraat  &  R.  E.  Setterington,  eds.  Publications 
Dept. ,  Electronic  Components  and  Materials  Division, 
Saugerties,  NY,  Sec.  Ed.,  Ch.  2,  4-22  (1974). 

3.  Martin,  G.  E.  "Vibrations  of  Coaxially  Segmented  Longi¬ 

tudinally  Polarized  Ferroelectric  Tubes,"  J.  Acoust.  Soc. 
Am.,  36,  1496-1506  (1964). 

4.  Martin,  G.  E.  "On  the  Theory  of  Segmented  Electromechanical 

Systems,"  J.  Acoust.  Soc.  Am.,  ^6,  1366-1370  (1964). 

5.  Cady,  W.  G.  Piezoelectricity.  McGraw-Hill,  NY,  Vol.  1,  Ch.  XIV, 

333-350  (1946). 

6.  Woollett,  R.  S.  "Transducer  Comparison  Methods  Based  on  the 

Electromechanical  Coupling  Coefficient  Concept,"  National 
Convention  Record,  I.R.E.,  _5,  23-27  (1957). 

7.  Cerpak,  V.  A.  "Dynamic  Lumped  Parameters  of  Composite  Piezo¬ 

electric  Transducers,"  Sov.  Phys.  Acoust.  ,  23,  //3,  246-250 
(1977). 

8.  Brickman,  A.  D.  "Preloading  of  Acoustic  Transducers  for 

High-Pressure  Operation,"  J.  Acoust.  Soc.  Am.,  34,  305-311 
(1972). 

9.  Tucker,  D.  G. ,  &  Gazey,  B.  K.  Applied  Underwater  Acoustics. 

Pergamon  Press,  Oxford,  154  (1966). 

10.  Glennite  Piezoceramics  Catalog.  Gulton  Industries,  Inc., 

Fullerton,  CA. 

11.  Ignizio,  J.  P.  Personal  Communication. 


AD-A075  282  PENNSYLVANIA  STATE  UNIV  UNIVERSITY  PARK  APPLIED  RESE--ETC  F/G  9/5 
THE  DESI6N  OF  PIEZOELECTRIC  TRANSDUCERS  USING  GOAL  PROGRAMMING. <U> 
JUN  79  D  F  MCCAMMO  N00024-79-C-6043 

UNCLASSIFIED  ARL/PSU/TM-79-140  NL 


DISTRIBUTION 


Commander  (NSEA  09G32) 

Naval  Sea  Systems  Command 
Department  of  the  Navy 

Washington,  DC  20362  Copies  1  and  2 

Commander  (NSEA  0342) 

Naval  Sea  Systems  Command 
Department  of  the  Navy 

Washington,  DC  20362  Copies  3  and  4 

Defense  Documentation  Center 
5010  Duke  Street 
Cameron  Station 

Alexandria,  VA  22314  Copies  5  through  16 


l 


l 


